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Simulating 

how  archaeological  observables  are  like  they  are 


INPUT  1 
INPUT  2 
INPUT  a 
INPUT  4 
INPUT  & 
INPU  T 6 
INPUT  7 


OUTPUT  1 
OUTPUT  2 
OUTPUT* 


We  build  a “model”  of  the  causal 
mechanism  associating 

human  activity-physics-observable  properties 

Inducction-Archaeological  Experimentation 


Facts 


Explaining 

why  archaeological  observables  are  where  they  are? 
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Simulating 

how  archaeological  observables  are  where  they  are 


We  build  a “model”  of  the  causal 
mechanism  associating 

human  activity-accumulation-spatiotemporal 

variation 

Inducction-Statistical  Modeling-Prediction 


The  coordinate 
chart  ctiK’  — * M 


The  manifold  M 
in  the  ease  that  k = 1 


Explaining 

why  human  society  changed  yesterday, 

changes  now  and 
it  will  change  tomorrow 


• The  change  from  Animality  to  Humanity 

• The  origins  of  Cultural  Diversity  and  Ethnogenesis 

• The  Origins  of  Social  Inequality 

• The  Origins  of  Agriculture 

• The  origins  of  the  State 

• The  Origins  of  Capital  Accumulation 
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Historians  explain  past  social  events  by 
showing  how  their  results  and 
consequences  fit  into  a causal  structure, 
that  is  to  say,  a vast  network  of 
interacting  actions  and  entities,  where  a 
change  in  a property  of  an  entity 
dialectically  produces  a change  in  a 
property  of  another  entity 
(transformation). 


annot  see  in  the  present  what 
ned  in  the  past.  We  would  need  a 
time  machine 


Explaining 

why  human  society  changed  yesterday, 

changes  now  and 
it  will  change  tomorrow 
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Bottom  line:  pay  attention  in 
tkifi  & wow  won’t  te  a£ 


There  is  no  way  of  actually  going  back 
into  the  past  to  test  a historical 
hypothesis. 


Explaining 

why  human  society  changed  yesterday, 

changes  now  and 
it  will  change  tomorrow 


If  we  cannot  enter  into  the  past  that  once 
existed,  we  can  enter  into  a virtual  past 
that  we  have  created  artificially  according 
to  what  we  know  for  certain  about  such 
past 


Simulating 

how  human  society  changed  yesterday, 

changes  now  and 
it  will  change  tomorrow 


Historians  seem  to  be  placed  in  a virtual 
world  extracted  from  a narrative, 
assumed  true,  written  by  an  individual 
who  saw  someone  doing 
something  in  the  past 


WBD  2010  illustrations  O Joe  Berger 


Simulating 

how  human  society  changed  yesterday, 

changes  now  and 
it  will  change  tomorrow 


Everyone 

SHOULD  BELIEVE 
IN  SOMETHING.... 


I BELIE  VE  I’LL  WRITE 
ANOTHER^CHA  PTER^ 


Such  virtual  worlds  are  expressed 
narratively,  using  verbal  language 


Simulating 

how  human  society  changed  yesterday, 

changes  now  and 
it  will  change  tomorrow 


isssMoiK  DrtniGF  O 

"IT'S  PEAL.  HERE  IS  A LETTER  OF 


AUTHENTICITY." 


In  archaeology  we  do  not  have  past 
witnesses,  nor  any  narrative  to  build  our 
own  “virtual”  worlds. 


Is  Archaeology  an  “impossible”  discipline 
beyond  what  ancient  materials  can 
suggest  from  their  material  properties? 


Simulating 

how  human  society  changed  yesterday, 

changes  now  and 
it  will  change  tomorrow 


It  has  been  suggested  that  Human 
Nature  has  not  changed  from  the  past  to 
the  present. 

Is  it  then  possible  to  understand  the  past 
from  what  we  know  human  behavior  and 
mentality? 

This  perspective  is  full  of  unprovable 
aphorisms 


Simulating 


how  human  society  changed  yesterday, 

changes  now  and 
'll  change  tomorrow 

Learning  by  Doing: 

Can  we  build  an  “artificial”  past? 

Let  us  consider  this  possibility: 


Simulating 

how  human  society  changed  yesterday, 

changes  now  and 
it  will  change  tomorrow 


Take  1000,  better  2000  humans. 

The  same  quantity  of  men  and  women. 
The  same  quantity  of  tall  and  short, 
clever  and  dull 


Simulating 

how  human  society  changed  yesterday, 

changes  now  and 
it  will  change  tomorrow 


Take  them  all  modern  technology 
(specially  smartphones!)  and  let  them 
build  a prehistoric  village,  and  reproduce 
the  same  activities  in  the  same  way  as 
our  ancestors 


Simulating 

how  human  society  changed  yesterday, 

changes  now  and 
it  will  change  tomorrow 


"You've  successfully  completed  the  transaction 
to  sell  your  soul.  A digital  receipt  will  be 
e-mailed  to  you." 


And  now  the  most  difficult  thing:  all  of 
them  must  sign  a contract  so  their 
children  only  learn  the  current  technology 
at  the  village  and  nothing  from  our 
present. 

And  repeat  this  contract  for  the  children 
of  their  children  for  the  next  4000 
years,  in  order  to  observe  changes  and 
transformations. 


Simulating 

how  human  society  changed  yesterday, 

changes  now  and 
it  will  change  tomorrow 


Is  this  as  “impossible”  as  it  seems? 

Instead  of  people,  let  us  program  “virtual 
People” 

Instead  of  a prehistoric  village,  let  us 
build  a virtual  village 

Instead  of  human  labour  and  behavior, 
let  us  write  some  algorithmic  functions 


Simulating 

how  human  society  changed  yesterday, 

changes  now  and 
it  will  change  tomorrow 


And  let  “virtual  agents”  live  in  a virtual 
environment  and  observe  how  everything 
evolves  and  changes  as  time  pass 


“Virtual”  means 


A reference  to  a “model”, 

a replication 

The  idea  of  something  acting  as  a 
surrogate 


The  idea  is  to  be  able  to“understand”  the  past 
because  we  have  “recreated”  it  inside  the  computer 
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What  is  a model? 


A representation  of 
some  (not  necessary  all) 
the  elements  of  an  entity 

The  aim  of  the  model  is 
to  allow  the 
understanding  of  the 
entity’s  structure  or 
behavior. 

It  allows  to  “experiment” 
with  a surrogate  of  that 
entity,  modifying  it 
according  to  the 
intended  final  goals 
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Simulating  Human  Society 


According  to  this  perspective,  the 
building  blocks  for  explanations  in 
the  social  domain  are 
products  (people,  goods,  informatior 

production  (human  labor,  social 
action),  and 

events  (the  context  in  which 
production  took  place) 

organized  such  that  they  are 
productive  of  some  changes  - 
regular  or  irregular-  from  start-up 

to  finish  or  termination  conditions 


Agent  Based  Modeling. 


Everything  is  an  “agent 


Environment 
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Memory 
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Resources 

• 

Decision  making  sophistication 

• 

Rules  to  modlify  behavioral  rules 

J 

Figure  1:  A typical  agent 


Agents  ->  Things  / Tools 


An  agent  simulating  a thing  is  just 
some  computer  code  with  specific 
attributes: 

Size 

Shape 

Texture 

Substance 

Mechanical  properties 


Agents  ->  Animals 


An  agent  simulating  an  animal  is  just 
some  computer  code  with  specific 
attributes: 

Size 

Sex 

Taxon 

Ethology 

Biomass 


Agents 


> Persons 


An  agent  simulating  a person  is  just 
some  computer  code  with  specific 
attributes: 

Age 

Sex 

Energy  needs 

Acquired  Energy 

Culture 

Knowledge 

Ethics  Values 

Goals 

Plans 


Agents->  Persons 


HOW 


TO  BUILD 
A PERSON 

A Prolegomenon 


The  idea  of  reproducing  people 
inside  a computer  is  not  very 
popular,  nor  it  is  well  accepted,  even 
at  scientific  domains. 

Do  you  play  with  people? 

Are  you  trying  to  predict  the  future? 

We  are  real  people , and  not  robots! 


An  Agent  Architecture 


26/04/2016 


31 


AG  E NTS 


REACTIVE 

AGENTS 


DELIBERATIVE 

AGENTS 


SIMULATING  HUMAN  DECISIONS 


Rational  Choice 


, . ' "Kill 

Jncertain 
World 


The  Psychology  of  Judgement 
and  Decision  Making 


Optimization 
Game  Theory 
Probability  Theory 
Fuzzy  Logic 
Heuristics 


SIMULATING  HUMAN  DECISIONS 


Optimization 

To  be  rational,  any  decision 
should  imply  a way  to 
evaluate  the  advantages  or 
drawbacks  of  each 
decision.  This  is 
implemented  in  terms  of 
utilities  to  be  maximized. 


minimize  /q  (ir) 

X 

subject  to  fi(x)  < h,  i = 1, . . . , m. 


SIMULATING  HUMAN  DECISIONS 


Player  2 
adopts 
strategy  A 

Player  2 
adopts 
strategy  B 

Player  1 

adopts 
strategy  A 

414 

1/3 

Player  1 
adopts 
strategy  B 

3/1 

2/2 

A sample  coordination  game  showing  relative 
payoff  for  player  1 / player  2 with  each  combination 


Optimization 

It  is  implemented  in  form 
of  Game  Theory. 

In  game  theory,  the  Nash  equilibrium  is 
a solution  concept  of  a 
non-cooperative  game  involving  two  or 
more  players,  in  which  each  player  is 
assumed  to  know  the  equilibrium 
strategies  of  the  other  players,  and  no 
player  has  anything  to  gain  by  changing 
only  their  own  strategy.  If  each  player 
has  chosen  a strategy  and  no  player  can 
benefit  by  changing  strategies  while  the 
other  players  keep  theirs  unchanged, 
then  the  current  set  of  strategy  choices 
and  the  corresponding  payoffs 
constitutes  a Nash  equilibrium 


SIMULATING  HUMAN  DECISIONS 


Institution  which  provides  frames  of  reference  (FR)  to 

the  gaps 


Information 
flow  filling 
the  gaps 


Information 
flow  from 
the  context 


Emerging 

state 


Agent -Agent  Relationships  (AAR) 
Agent-Environment  Relationships  (AER) 

Environment  Laws  (EL) 


Planned 
state 


Optimization 

To  be  rational  and  be 
capabe  of  maximizing  or 
minimizing  utilities,  an 
agent  should  learn  the 
evolving  environment  of 
landscape  features  and 
actions  performed  by  other 
agents  in  the  neighborhood 


Environment 


SIMULATING  HUMAN  DECISIONS 


COMPLEX  PROBLEMS 


COGNITIVE 

LIMITATIONS 


INFORMATION 

IMPERFECTION 


TIME 

CONSTRAINTS 
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Incomplete  Information 


SIMULATING  HUMAN  DECISIONS 


Interpreting  the  Fuzzy  inference  Diagram 


Fuzzy  Logic  to  implement 
Bounded  rationality  mechanisms 


SIMULATING  HUMAN  DECISIONS 


RES  QURCE  AVAILAB  LITY 
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Beyond  the  individual  level,  accumulated  decisions 
can  be  explored  probabilistically 


And  Agent  Based  Models  are 
systems  ef  connexted  expert 

systems 


Agents  ->  Sets  of  Persons 


Because  in  Archaeology  the 
Individual  agent  is  usually  invisible,  it 
Is  very  common  to  program  an  “artificial 
Society”  at  the  level  of 
HOUSEHOLD  (family) 

LOCAL  GROUP  (Settlement) 

NATION 

Attribute  set  is  then  different. 


Virtual  Agents  in  a Virtual  Environment 


This  is  the  easiest  part, 
Because  we  can  integrate  GIS 
Models  into  an  ABM. 

You  can  use  a vectorial 
Geographical  database 
Or  a raster-based  GIS 
(more  usual) 


Virtual  Agents  in  a Virtual  Environment 


GIS  environmental  data 

Temperature 
Crops 

Precipitation 


Cellular  automaton 

Pest  survival  (T  ; P ; crops) 
Pest  dispersal  (density) 
Pest  fecundity  (T) 


GIS  social  data 

Village  location 
Population  size 


Agents 

Farmers  - pest  control  knowledge 

- team 

- infestation  state 


Virtual  Environments  are  usually  divided  into  “patches”,  that  is  homogenous 
Areas  with  the  same  quantity  of  resources,  vegetation,  or  geographical  features 


Virtual  Agents  in  a Virtual  Environment 


Alison  J.  Heppenstall 
Andrew!  Crooks 

Linda  M.  See  * Michael  Batty  Editors 

Agent-Based 
Models  of 
Geographical 
Systems 


Springer 


It  is  up  to  you  to  decide  to 
Program  the  Environment 
As  a particular  entity  with  values  (GIS) 
Or  as  a series  of  agents: 

For  instance: 

GIS:  Number  of  Deer  (Attribute  Patch) 
Deers  are  Agents  within  the  system 


Agent  Based  Dynamics 


► An  agent  is  ‘alive’  in  its  environment.  An  agent 
has  direct  interactions  with  its  world. 

► The  agent  may  act  on  the  environment,  which  in 
turn  provides  perceptions  to  the  agent. 

► The  complexity  of  these  interactions  is  primarily 
constructed  in  the  agent,  and  not  in  the 

environment.  An  agent  observes  and  interprets 
the  world.  The  environment  may  change,  but  the 
agent  will  have  to  observe  these  changes. 


Agent  Based  Dynamics->  Environment 


Plant  biogeography  and  vegetation  dynamics 
(resource  competition,  seasonality  and  growth) 


Fundamental  niche 


Species  A, 


A,  dominates/ 

B,  coexists  « 

B,  dominates 

\Vn<*  C,  coexists  c dominates 
I B.  coexists 

Environmental  space  (1) 


& Environmental  space  (n) 


a.  Environmental  space  (3) 


a Environmental  space  (2) 


Plant  and  soil  biophysics  and  biogeochemistry 
(carbon,  water  and  nutrient  budgets) 


H,0 


You  should  introduce  “natural”  dynamics 
clima 
pluviosity 

natural  reproduction 
Natural  catastrophes 


To  build  the  “environmental  engine”, 
you  use  paleoecological  data  and 
ecological  mechanisms 


Nutrients  " 


Agent  Based  Dynamics->  People 


COHTHiilirPONi  TO  ECONOMICS 


I'mni  Msm  pilhtn 
PFEkJMPtZ 


Agent-Based 

Computational 


Ki  Tunprrtry^  C^i-r  Urderatindmfl: 
uf  PoiHHgraptiG  Bebavigur 


R7 

It-.!.-  IMfrltrifEMD 


DEMOGRAPHY: 

Demographical  mechanisms  are  relatively 
Well  known. 

We  should  program  how  agents  die 
and  new  agents  born,  and  how  the 
number  of  living  agents  evolves. 

This  is  not  just  biological  reproduction. 

We  can  introduce  social  norms  to 
constrain  human  reproduction 


Agent  Based  Dynamics->  People 
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Thinking 
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Acting 
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Logical  Lo  Linda  tin  ns  tor 

o 

Rational  Decision  Making 


Behavioral  Mathematics 
for  Game  Al 


The  Mathematics 
of  Behavior 


Agents  behave: 
move 
hunt 

gather  vegetables 
produce  subsistence 
share 
interact 
kill 
die 


Earl  Hunt 


Coding  "Human  action" 

to  move 

ask  families  [ 

set  - neighborhood - and - group 
move -to  one -of  my -neighborhood 
set  heading  random  360  fd  0.5  ] 

;if  movement  - rule  = "follow  local  leader" 

; [ ask  families  [ 

; set -neighborhood -and -group 

; set  leader  one-of  (turtle- set  self  my-group)  with-max 
[capacity] 

; move -to  one-of  [my -neighborhood]  of  leader  ] 

;] 

End 


What  is  a “Social  Mechanism 


Mechanisms 

are  entities  and 
activities  organized 
such  that  they  are 
productive  of  regular 
changes  from  start  or 
set-up  to  finish  or 
termination  conditions 


Rwng-faChw 
Rrt*rtjiL.Milhlnn  £tfws 

Mechanism  and 
Causality  in  Biology 
and  Economics 


What  is  a “Social  Mechanism” 


Obviously,  the  word 
“mechanism”  is  here  a 
parable  of  how  social 
intentions,  goals  and 
behaviors  are  causally 
connected. 

It  should  explain  how 
social  activity  worked, 
rather  than  why  the  traits 
contributing  to  these 
activities  or  workings  are 
there 


Organizational 

Influences 


Latent 


Failures 


Unsafe 

Supervision 


Latent 


Failures 


4 


Preconditions 
for  Unsafe 
& Acts 


Latent 


Failures 


Unsafe 

Acts 


Failed  or 
Absent  Defenses 


Impact  of  Error 


What  is  a “Social  Mechanism” 


Is  there  any  human  aspect  that  cannot 
can  be  “simulated”  mechanically? 
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SOCIAL  INTERACTION 


Towards  a theory  of  “action” 


A Theory  of  Social  Action  (Leont’iev  1974) 

Social  action  can  be  defined  in  terms  of  purposeful 
changing  of  natural  and  social  reality.  Social  actions 
are  goal-directed  processes  that  must  be  undertaken 
to  fulfill  some  need  or  motivation.  They  are  conscious 
(because  one  holds  a goal  in  mind),  and  different  actions 
may  be  undertaken  to  meet  the  same  goal. 

However,  an  action  can  be  an  intentional  action  without 
the  actor  having  to  be  aware  of  the  intention  from  moment 
to  moment.  Motivations  or  intentions  are  not  just  conditions 
for  developing  cognitive  activity,  but  real  factors  influencing 
productivity  and  perceivable  structure. 


Theory  of  Agency  (Bourdieu,  Giddens) 

"agency"  refers  to  the  capacity  of 
individual  humans  to  act  independently 
and  to  make  their  own  free  choices. 

Reaction  against  the  idea  that  our  social  existence 
is  largely  determined  by  the 
overall  structure  of  society. 

This  theory  sees  structure  and  agency  as 
Complementary  forces  - structure  influences 
human  behaviour, 

and  humans  are  capable  of  changing  the 
social  structures  they  inhabit. 


SOCIAL  INTERACTION 


Towards  a theory  of  “action” 


A Theory  of  Social  Action  (Leont’iev  1974) 

Social  actions  cannot  be  understood  without  a frame  of  reference  created  by  the  corresponding  social  motivation 
or  intention.  Leont’ev  describes  social  activity  as  being  composed  of  subjects,  needs,  motivations,  goals, 
actions  and  operations  (or  behavior),  together  with  mediating  artifacts  (signs,  tools,  rules,  community,  and 
division  of  labor) 

• A subject  is  a person  or  group  engaged  in  an  activity. 

• An  intention  or  motivation  is  held  by  the  subject  and  it  explains  activity. 

• Activities  are  realized  as  individual  or  cooperative  actions.  Chains  and  networks  of  such  actions  are  related  to 

each  other  by  the  same  overall  object  and  motivation. 

• Actions  consist  of  chains  of  operations,  which  are  well  defined  behaviors  used  as  answers  to  conditions  faced 

during  the  performing  of  an  action. 

• Activities  are  oriented  to  motivations,  that  is,  the  reasons  that  are  impelling  by  themselves. 

• Each  motivation  is  an  object,  material  or  ideal,  that  satisfies  a need. 

• Actions  are  the  processes  functionally  subordinated  to  activities;  they  are  directed  at  specific  conscious  goals. 

• Actions  are  realized  through  operations  that  are  the  result  of  knowledge  or  skill,  and  depend  on  the  conditions 

under  which  the  action  is  being  carried  out. 


Social  Influence  Theory 

M Lewenstein,  A Nowak,  B Latane  - Statistical  mechanics  of  social  impact  Physical  Review  A,  1992 


I = f(  SIN). 
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Intention 
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Social  impact  appears  to  be  the  result  of  a multiplicative  function  based 
on  the  number  of  people  acting  on  the  target.  Within  this  equation,  the 
strength  (S)  is  a measure  of  how  much  influence,  power,  or  intensity  the 
target  perceives  the  source  to  possess. 

Immediacy  (/)  takes  into  account  how  recent  the  event  occurred  and 
whether  or  not  there  were  other  intervening  factors. 


Behavior 


Actual 

Behavioral 

Control 


Furthermore, 

/ = sNp 

That  is,  some  power  (p)  of  the  number  of  people  (A/)  multiplied  by  the 
scaling  constant  (s)  determines  social  impact. 


COMPUTER  SIMULATION 


Why  do  I what  I do  the  way  I do? 


Because  other  people  do 
the  things  they  do 
in  the  way  they  do 


COMPUTER  SIMULATION 
Why  am  doing  this? 

Because  other  people  do  the  things  they  do  in  the  way  they  do” 


I’M  NOT  DOING  WHAT  MY  NEIGHBOUR  DOES 

HOWEVER,  MY  RATIONAL  AND  IRRATIONAL  DECISIONS 
ARE  AFFECTED  BY  THE  WAY  MY  NEIGHBOURS 
TAKE  RATIONAL  AND  IRRATIONAL  DECISIONS 


SOCIAL  INTERACTION 


©Original  Artist 

Reproduction  rights  obtainable  from 
www  C artoonStock  .com  — , 


SOCIAL  INTERACTION 
Towards  a theory  of  “action” 


Theory  of  Agency  ( Beyond  Bourdieu,  Giddens) 
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RUNING  THE  SIMULATION 


Let  us  execute  the 
Model: 


Agents  “act”, 

They  move, 

They  use  energy 
They  acquire  energy 
They  reproduce 
They  interact  with  others 
They  change  their 
previous  knowledge 


RUNING  THE  SIMULATION 


To  implement  activity  we  have  to 
specify  the  order  in  which  events 
happen. 

A cycle:  all  an  agent  shoul  do  in 
one  day, 
one  season 
one  year 

one  hundred  years 


Predator-Prey  classical  model 


RUNING  THE  SIMULATION: 


Simple  State 

Control  is  always  located  in  one  of  the  simple 
states.  States  have  Entry  and  Exit  actions 

Transition 

Can  be  triggered  by  an  external  or  internal 
event,  condition  or  timeout.  Has  Action 

Composite  State 

A group  of  states  with  common  behavior 
History  Pseudo  State 

Denotes  last  visited  state  in  the  composite  state 
Branch  Pseudo  State 

Specifies  conditional  branching  of  transitions 


RUNING  THE  SIMULATION:  The  Idea  of  Time 


As  soon  as  attribute  and  global 
Variables  change  their  values 
We  obtain  an  idea  of  time. 

Time  as  the  sequence  of  changes 


RUNING  THE  SIMULATION:  emergent  behavior 


But  “changes”  are  not 
easy  to  detect,  because 
They  are  not  sequential 
ror  linearly  related: 

EMERGENCE  is  a process 
whereby  larger  entities, 
patterns,  and  regularities 
arise  through  interactions 
among  smaller  or  simpler 
entities  that  themselves  do 
not  exhibit  such  properties 
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Complex  Adaptive  Behavior 
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Simple  Self-Organized 
Local  Relationships 
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MONITORING  THE  SIMULATION 


Edited  by  Pierre  Demeulenaere 


CopyrluMtd  Maori* 


The  Oxford  Handbook  of 

ANALYTICAL 

SOCIOLOGY 


Analytical  Sociology 
and  Social  Mechanisms 
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Experiments 
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We  should  create  numeric  output  variables  to 
“measure” 

the  system  dynamics. 

THE  IDEA  OF  “SOCIAL  INDEX” 

It  is  a composite  measure,  i.e.  measurements  based 
on  multiple  data  items. 

Three  common  composite  measures  include: 
indexes  - measures  that  summarize  and  rank  specific 
observations  (usually  on  the  ordinal  scale). 
scales  - advanced  indexes  whose  observations  are 
further  transformed  (scaled)  due  to  their  logical  or 
empirical  relationships; 
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SimulPast 

....simulating  the  past  to  understand  human  behaviour 
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Archaeologists 
have  been  doing 
material 
classifications  to 
discover  "cultural" 
limits  and 
differentiation,  as 
surrogates  of 
different  ethnic 
groups,  avoiding  a 
careful  analysis  of  the 
very  meaning  of 
“ethnicity” 


Regions: 

1.  Eastern 

2.  Central 

3.  Western 

4.  Iberian 

5.  Italian 
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SOCIAL  VARIABILITY 


Design  Concepts: 

Ethnicity  as  Cultural  Differentiation 


An  ethnic  group  is  a group  of  people  whose  members 
explicitly  regard  themselves  and  are  regarded  by  others 
as  truly  distinctive,  through  a common  heritage  that  is 
real  or  assumed-  sharing  “cultural”  characteristics. 

The  emergence  of  such  identification  is  the 
consequence  of  different  historical  processes  and  social 
mechanisms,  which  take  place  within  the  webs  of  group 

and/or  individual  interaction. 


Design  Concepts: 

Ethnicity  as  Cultural  Differentiation 


The  question  is  “why  groups  of  people  are  the  way  they 
are”  in  terms  of  how  they  acted  within  a social  aggregate 
their  previous  activity  contributed  to  build.  The  complex 
interplay  of  social  actions,  people  and  the  consequences 
of  their  actions  explain  ethnicity  by  showing  how  social 
aggregation  fit  into  a causal  structure,  that  is  to  say,  a 
vast  network  of  interacting  actions  and  entities,  where  a 
change  in  a property  of  an  entity  dialectically  produces  a 
change  in  a property  of  another  entity 


Design  Concepts: 

Ethnicity  as  Cultural  Differentiation 


Cultural  Standartizacion  as  an  Emerging 
Phenomenon 

As  a consequence  of  different  forms  of  social 
interaction  and  exchange  between  related 
individuals  cultural  consensus  may  emerge, 
and  communities  of  people  can  aggregate  into 
greater  groups  culturally  homogenous,  what 
affect  social  reproduction,  increasing  similarity 
in  the  long  run,  and  reducing  the  risk  of  being 
attacked  by  groups  identified  as  "enemies" 

(out  of  the  new  cultural  consensus). 


Design  Concepts: 

Ethnicity  as  Cultural  Differentiation 


ethnogenesis  (at  the  macro-scale)  appears  as  a 
consequence  of  a long  period  of  repeated  constraints  to 
cooperation  due  to  stable  identity  formation  (at  the 
micro-scale) . 

Both  are  a contradiction  between  social  inertia 
(knowledge  inheritance)  and  cultural  consensus  (social 
similarity)  built  during  cooperation  and  labor  exchange. 


Design  Concepts: 

Ethnicity  as  Cultural  Differentiation 


Schelling  proposed  a model  of  a 
population  randomly  displaced  on  a two 
dimensional  grid.  The  grid  hosts  the 
agents  and  each  agent  has  a neighbor 
made  up  of  four  cells,  those  immediately 
up,  down,  left  and  right  to  the  agent. 

Each  agent  is  happy"  or  not  depending 
on  the  number  of  its  neighbors:  if  the 
percentage  of  neighbors  of  the  same 
breed  over  the  total  number  of  neighbors 
is  under  a fixed  threshold,  the  person  is 
unhappy.  Agents  move  in  order  to 
enhance  the  probability  of  being  happy 
with  most  neighbours.  Schelling  shows 
that  after  a limited  number  of  moves,  a 
segregation  pattern  emerges,  and  that 
this  macro"behavior  is  largely 
unpredictable  based  only  on  the 
only  simple  rule  each  agent  follows, 
namely  that  to  move  if  it's  unhappy. 


SCHELLING  Segregation 


powered  by  NetLogo 


Design  Concepts: 

Ethnicity  as  Cultural  Differentiation 


AXELROD’S  Diffusion  of  Culture 


Axelrod's  model  of  cultural  dissemination.  It 
models  a population  of  actors  that  hold  a 
number  of  cultural  attributes  (called  features) 
and  interact  with  their  neighbors.  Dynamics 
are  based  on  two  main  mechanisms.  First, 
agents  tend  to  chose  culturally  similar 
neighbors  as  interaction  partners  (homophily). 
Second,  during  interaction  agents  influence 
each  other  in  a way  that  they  become  more 
similar.  The  interplay  of  these  mechanisms 
either  leads  to  cultural  homogeneity  (all  agents 
are  perfectly  similar)  or  the  development  of 
culturally  distinct  regions.  The  model  allows 
studying  to  which  degree  the  likelihood  of 
these  two  outcomes  depends  on  the  size  of 
the  population,  the  number  of  features  the 
agents  hold,  the  number  of  traits  (values)  each 
feature  can  adopt  and  the  neighborhood  size 
(interaction  range). 


Design  Concepts: 

Ethnicity  as  Cultural  Differentiation 


DIFFERENCES  WITH  OUR  MODEL: 

1)  Integrating  Axelrod  and  Schelling:  mobility  is  allowed,  and  transmission  of 
cultural  features 

2)  Culture  is  not  fixed,  but  evolve  as  soon  as  cooperation  appears  to  be 
advantageous 

3)  Random  cultural  drift  simulates  internal  process  in  addition  to  external 
interaction 

4)  Movement  is  not  random,  but  depends  on  survival 

5)  Survival  depends  on  cooperation 

6)  Cultural  transmission  depends  also  on  survival  and  hence  on  cooperation 

7)  Cooperation  and  reciprocity  depends  on  the  probabilities  of  survival 


Purpose: 

how  diversity  and  self-identification  can  emerge  in  small-scale 

societies 


In  this  preliminary  and  simplified  computer  simulation  we  are 
exploring  the  consequences  that  labor  exchange  and  territorial 
mobility  in  an  artificial  unconstrained  world  has  on  identity 

formation  and  negotiation. 

We  expect  to  be  able  to  discern  if  cultural  diversity  emerge  as  a 
result  of  social  decisions  only,  or  if  it  is  the  result  of  constraints  on 
mobility  generated  by  geography  and  the  irregular  distribution  of 

resources,  both  in  space  and  time. 


Purpose: 

how  diversity  and  self-identification  can  emerge  in  small-scale 

societies 


Agents  should  be  able  to  update  and  consequently  modify  their  actual  identity 
(“culture”)  when  what  they  inherited  proves  to  be  not  useful  to  cooperate. 

What  agents  have  inherited  at  birth  is  usually  a central  tendency  of  parents 
identity.  Given  that  the  process  of  mating  and  selecting  reproductive  partner  is 
socially  and  politically  mediated,  the  cultural  consensus  to  be  transmitted 
depends  on  the  way  such  cultural  consensus  is  built  at  the  level  of  the 
reproductive  unit,  and  hence  on  the  social  and  political  origins  of  the 
reproductive  unit  itself  (“family”). 
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SPATIAL  INDEX 


Social  fractionalizahon  (ELF) 

Depth  of  differences 

Generalized  fractionalizahon 

Cultural  polarization  (RQ) 

spatial-index 
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Number  of  family  deaths  (current  tick) 
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Mean  labor  per  family 
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Mean  technology  of  families 
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COOPERATIVE  UNITS  STATISTICS 


display-resorces? 

initial-population  150 

mean-resource-on-patc...  7000 

average-technology  1.05 

labor-average  6 

movement  35 

■ ■ 

print-spatial-matrix? 

sd-resource-on-patches  1000 

diversity  0.1 

average-storing-factor  0.5 

internal-change-rate  0.01 

Number  of  agents  in  cooperative  units 

Number  of  isolated  agents 

139 
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Number  of  cooperative  units 

Number  of  agents  in  the  largest  cooperative  unit 
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Avg.  similarity  threshold  in  cooperative  units 

Std.  dev.  of  similarity  thres.  in  cooperative  units 

0.0030 

0.02 

ENERGY  $ 

3TATISTICS 

Std.  deviation  of  max-resource 

1008.552934914555 

Total  collected  energy 

Standard  deviation  of  collected  energy 

237804.55 

555.48 

364.17 

What  is  possible  within  the  model? 
Exploring  the  Parametric  Space 


What  is  possible  within  the  model? 

Exploring  the  Parametric  Space 
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Number  of  agents  in  cooperative  units 

Number  of  isolated  agents 

139 

7 

Number  of  cooperative  units 

Number  of  agents  in  the  largest  cooperative  unit 

5 

104 

Avg.  similarity  threshold  in  cooperative  units 

0.0030 

Std.  dev.  of  similarity  thres.  in  cooperative  units 

0.02 

ENERGY  5 

STATISTICS 

Std,  deviation  of  max-resource 

1008.5529349H555 

Total  collected  energy 

Standard  deviation  of  collected  energy 

237804,55 

555.48 

Collected  energy  per  labor 

364.17 

The  Model 


Do  the  degree  of  cooperation  and  its  social  consequences  (identity  variation  at 
the  short  term)  affect  global  identity  variation  at  the  long  term? 


The  model  has  been  Created  to  predict  degree  of  Cultural  Diversity  and 
Ethnogenesis,  which  is  measured  according  five  indexes: 


1.  Social  Fractionalization 
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2.  Depth  of  Differences 

Non-normalized  Euclidean  distances  between  identity  vectors. 


3.  Generalization  of  Fractionalization 
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4.  Generalized  Resemblance 
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5.  Generalized  Spatial  Aggregation  and  Territoriality 
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Entities,  state  variables  and  scales: 

Agents 


agents  are  modeled  as  a number  of 
individuals  acting  as  a decision  unit, 
that  is,  what  we  usually  name  “family”  in 
real  live,  and  what  anthropologists 
define  as  “domestic  unit”  or  “household” 

Agents  are  obliged  to  constant 
geographical  mobility,  given  their 
extreme  dependence  to  local  carrying 
capacity  and  diminishing  returns  from 
labor  given  that  they  cannot  restore 
what  they  have  extracted  from  nature 


Entities,  state  variables  and  scales: 

Agents 


AGE  (number  of  lived  time-steps  at  each  time-step) 
MAX-AGE.  A condition,  after  what  the  agent  dies. 
Fixed  for  all  agents 

SURVIVAL.  Probability  of  survival  for  each  agent.  It 
depends  on  the  actual  value  of  energy.  It  determines 
if  the  household  survives  or  dies. 

LABOR  UNITS:  (/,)  (a  Poisson  distributed  parameter 
counting  the  aggregated  quantity  of  labor  from  all 
groups) 

TECHNOLOGY  (b):  A parameter  representing  the 
aggregated  efficiency  of  labor  obtained  when 
increasing  the  number  of  manufactured  tools).  It 
starts  at  1 (lack  of  tools)  and  has  not  an  upper 
maximum. 

ENERGY  = ( ei}  Produced  food,  expressed  in 
kilocalories. 

SURPLUS  (S|)  The  difference  between  energy 
produced  and  energy  consumed.  It  is  storted  for  later 
use. 

SURVIVAL  THRESHOLD,  (e,):  in  terms  of  the 
quantity  of  calories  all  agents  included  in  an  agent 
that  represents  a regional  group  of  local  groups  need 
to  be  able  to  live  a season  long  (six  months). 


Entities,  state  variables  and  scales: 

Agents 


CULTURE:  A 10x10  matrix  of  natural  number.  At 
the  beginning  of  the  simulation  all  agents  have 
the  same  culture.  It  is  a measure  of  IDENTITY 
MY-GROUP  A list  of  connected  households  in  the 
same  social  aggregate. 

MY-NEIGHBORHOOD  (list  of  agents  around) 
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Entities,  state  variables  and  scales: 

CULTURAL  IDENTITY 
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Entities,  state  variables  and  scales: 

CULTURAL  IDENTITY 
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Entities,  state  variables  and  scales: 

CULTURAL  IDENTITY 
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Organized  by  motivational 
similarities  and  dissimilarities 
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Entities,  state  variables  and  scales: 

Environment 


[ Interface  Information  Procedures 

O view  updates  | t v 

faster  ! on  ticks  I &ED 


Uniform  Distribution  of  Kilocalories: 

-between  100  and  15000  POOR  WORLD 

*between  100  and  20000 

-between  100  and  25000 

-between  100  and  40000 

-between  100  and  50000  RICH  WORLD 


We  have  built  a 
world  without 
topographical 
barriers,  where 
resources  are 
irregularly 
distributed 
across  geographical 
space,  and  with  a 
founding  population 
having  a single 
homogenous 
identity,  with  a 
constant,  but 
random  internal 
change  rate. 


Process  overview 


Let  us  imagine  a population  of  virtual 
agents,  moving  randomly  in  search 
for  resources,  and  organized  in 
households  . 

If  the  agent  survives,  and  there  are 
other  agents  in  the  neighborhood  with 
enough  cultural  consensus,  linked 
agents  have  the  opportunity  to  give 
birth  to  a new  agent  (reproduce)  and 
to  transfer  a combination  of  the 
identity  of  linked  agents  to  the  new 
generation. 


Process  overview 


Agents:  households 


Agents  all  have  a similar  identity  vector, 
and  they  begin  working  at  nearby  places. 
Very  soon  they  begin  to  move  to  avoid 
diminishing  returns  derived  from  local 
carrying  capacity.  When  moving,  they  lose 
contact  with  brothers  and  cousins,  and  their 
identity  began  to  change  in  isolation  given 
internal  mechanisms.  The  greatest  the 
geographical  distance,  the  less  the 
probability  of  renewing  contacts  and 
maintaining  cultural  consensus.  If  they 
randomly  find  a relative  in  their  peripatetic 
displacements,  a common  identity  is  rebuilt, 
but  it  is  a new  identity  and  not  the  previous 
one,  integrating  some  of  the  new  elements 
that  may  have  emerged  in  isolation. 


SURVIVING: 

Looking  for  food 


DYING: 

When  resources  are  not  enough 


MOVING: 

If  resources  are  not  enough 
the  agent  moves  to  a near  place 


When  cultural  consensus  is  low 
labor  is  performed  individually 


INDENTIFYING: 

Culturally  similar  people  within 
the  neighborhood 


COOPERATING: 

When  cultural  consensus  is  enough, 
labor  is  collectively  performed 


Survival  Cycle 


Inherent  identity 
at  birth 


Dyng: 

- If  resources  are  not  enough 
-If  they  reach  the  maximun  age 
-Randomly  30%  of  birth  are  discarded 
-Randomly  5%  of  agentes  are  descarded 


Evolution 
of  Identity 
Cycle 


Use  its  own  identity 
To  identify  the  others 


Every  10  cycles,  2 neighbour  agents 
Produce  a new  agent  with  an  identity 
That  combines  their  parents  identity 


If  cultural  consensus  is 
similarity  increase 


If  cultural  consensus  is  low, 
identity  mute  randomly 


high, 


Process  Overview.  Sequence  of 


Subprocesses.  SURVIVE 


SURVIVE:  Agents  start  with  enough  energy  to  survive  the  first  time-step. 
Energy  at  start-up  is  calculated  as 


"Eat  my  vegetables?  — I thought 
we  were  supposed  to  be  predators!1 


Subprocesses 


to  survive 

; harvest  and  hunting 
ask  families  [ 

let  group  - capacity  ( capacity  + sum  [ capacity  ] of  my- group  ) 

let  family - output  ( capacity  * ( group  - capacity  A ( returns  - to  - cooperation  - 

1 ) ) ) 

set  energy  ( ( energy  * (1  - energy - depreciation)  ) - subsistence  + family- 

output  ) 

set  energy  max  (list  0 energy  ) 

] 

end 


keGj(t) 


keGj(t) 


keGj(t) 


6-1 


with  0 > 1 


Subprocesses 


Each  agent  begins  by  using  existing 
surplus,  produced  at  the  end  of  the 
previous  time-step.  If  stored  surplus  (s,) 
is  equal  to  the  survival  threshold  (e(),  the 

agent  survives.  If  stored  surplus  is  not 
enough,  the  agent  should  produce  food 
(e,).  After  survival,  the  amount  of  food 

that  is  not  consumed  is  converted  into 
surplus. 


SURVIVE 


///'>,/  y\ 


Subprocesses.LOOK  FOR  FOOD 


Food  is  obtained  by  agent  /'  by  means  of  labor  (lt(t)) 
with  the  contribution  of  its  own  technology  b,( t), 
used  to  compensate  the  local  difficulty  of  producing 
food  (/?,).  In  this  simulation,  we  assume  that  any 
individual  can  produce  360  kilocalories  in  a single 
six-months  time-step  (2000  calories  per  day  * 180 
days). 


e,=  [360  l,(t)l  hi  ] * bi t) 

(/?,)  is  a Poisson  distributed  parameter  counting  the 
quality  of  soil  and  the  availability  of  water  and 
temperature:  the  poorer  quality  of  soils  and  the 
scarcer  is  water,  the  more  labor  or  more  technology 
is  needed  to  obtain  resources  up  to  survival 
threshold).  This  parameter  is  initiated  at  start  up  (a 
random  number  following  a Poisson  distribution 
whose  I is  a free  parameter  selected  by  the  user  at 
the  beginning  of  the  simulation),  and  changes  every 
time-step,  in  such  a way  that  at  odd  cycles  (warm 
season)  it  is  the  half  that  at  even  cycles  (cold 
season). 


* were  HUMTER5  AMD  GATHERERS  .BUT.  TOJKqHT 
WE  DECIDED  TO  EAT  OUT 


#20900081 


Subprocesses.HUNT-and-GATHER 

Energy  is  obtained  by  agent  / by  means  of  labor,  with  the  contribution  of  its  own  technology  {(5). 


1 + 


1 

1 


This  equation  measures  the  ability  to  obtain  resources  according  to  each  agent’s  individual  ability. 


Subprocesses. IDENTIFY- 


Many  procedures  depend  on  the  actual  degree  of  cultural  consensus  between  agents.  It  is 
calculated  in  terms  of  the  similarity  between  the  identity  of  agents.  The  identity  of  an  agent  is 
expressed  in  a binary  vector  including  the  presence/absence  information  of  cultural  elements, 
which  can  be  also  identified  in  the  archaeological  record.  Given  that  each  agent  is  in  fact  an 
integration  of  local  groups  (archaeological  sites:  settlements,  hillforts,  workshops, 
cemeteries,  ritual  places,  etc.)  such  a vector  contains  the  statistical  mode  of  all  attributes 
present  in  the  buffer  zone  corresponding  to  each  “regional”  agent. 


How  do  the  agents  interact?: 


• Two  agents  regard  themselves  similar  if  they  share,  at 
least,  a proportion  of  6 cultural  traits  (cultural  proximity) 


Let's  cooperate! 


Subprocesses.IDENTIFY- 

NEIGHR™tr<; 

Hamming  distance  = 3 


When  identifying  other  agents  in  the  area,  each  agent  calculates  a Hamming  distance 
between  each  pair  of  identity  vectors.  The  Hamming  distance  between  two  strings  of  equal 
length  is  the  number  of  positions  at  which  the  corresponding  symbols  are  different. 


} 


In  the  equation  dHAD  is  the  Hamming  distance  between  the  objects  / and  j,  k is  the  index  of  the 
respective  variable  reading  y out  of  the  total  number  of  variables  n. 


Subprocesses.  IDENTIFY 


How  big  should  be  this  threst^S^JlOT/ 
cooperation  between  agents  or  increase  the 
probabilities  of  being  attacked?  If  an  agent  is  able 
to  produce  enough  food  by  itself  and  survive,  it 
does  not  bother  to  identify  other  agents  in  the 
territory.  Otherwise,  agents  calculate  the 
percentage  of  consensus  needed,  depending  on 
how  much  they  need  food  or  tools  from  others  to 
survive.  The  greater  the  ease  with  which  an  agent 
obtains  needed  resources,  the  more  predisposed 
to  help  at  no  cost.  This 


st =— ([.  ,*l 

ICC  v 


A similarity  threshold  (ST)  expressing  the 
degree  of  independence  an  agent  has  on 
others 


Subprocesses : RATIONAL 


If  the  current  measure  of  i 
and  j is  above  the  calculated  value  of  ST(t)then  there  is  option 
to  decide  whether  exchange  or  not,  according  to  expected 
benefits. 

To  be  rational,  the  decision  whether  sending  the  requested 
food  or  refusing  the  proposed  exchange  should  imply  a way  to 
evaluate  the  advantages  or  drawbacks  of  this  behavior. 

Prisoner's  Dilemma 


Region  2 ( with  food  in  excess) 

Exchange 

Refuse  Exchanging 

Regionl 

(in  need) 

Exchange 

5,3 

0,5 

Refuse 

Exchanging 

5,0 

5,0 

In  our  case, 

Temptation  (T)  > Reward  (R)  > Punishment  (P)  = Sucker  (S) 

And  it  is  easy  to  see  that 

2R  > T + S 


Consequences  of  cooperation: 


Higher  output: 
— Individually: 


get  more  output  than  if  I act  individually: 

— Cooperation: 


0,(t)  = c, 


°i(t)  = ci(  Z Ck) 

\keGj(t)  / 


with  0 > 1 


O:  Increasing  returns-to-cooperation 


Subprocesses.  Cooperate 


Provided  there  is  enough  CULTURAL 
CONSENSUS,  agents  who  cannot  survive  on 
their  own  are  helped  by  “culturally  similar” 
neighbors  with  enough  amounts  of  labor  and 
technology 


1 


1 + 


1 


where  Gft)  is  the  total  amount  of  labor  the  group 
of  agents  that  cooperate  with  agent  / and  dp  ft) 

the  maximum  technology  within  the  group.  There 
is  an  additional  parameter  modifying  the  total 
effect  of  aggregated  labor  at  the  social 
aggregate  (Oft)),  capturing  the  idea  that 

cooperation  is  less  needed  when  there  are 
plenty  of  resources.  In  other  words,  it  measures 
the  added  value  that  cooperation  brings  to 
production  returns. 


MWe  are  neither  hunters  mr gatherers.  We  are  aaamtmts, ' 


Subprocesses.  UPDATE  "Culture" 


Once  the  agent  gets  enough  resource  for  its  own 
survival  (with  or  without  the  help  from  others),  the 
identity  vector  (CULTURE)  is  updated  towards  the 
statistical  mode  of  the  groups  identity. 

With  a fixed  probability  level  (95%)  each  agent  copies 
the  statistical  mode  of  identities  within  the  group. 

There  is  an  additional  source  of  identity  change, 
implemented  in  form  of  an  INTERNAL  CHANGE  RATE 
(IRC).  This  is  a random  value  (from  0 to  1,  usually  very 
small)  defined  in  analogy  to  the  probabilities  of 
internal  change  (invention,  mutation,  catastrophe, 
sudden  change). 


TECHNOLOGY  IS  ALSO  UPDATED  IN  THE  SAME  WAY, 
towards  the  most  efficient  value  in  the  emerging 
aggregate  of  cooperating  agents 


Subprocesses 


to  evolve  - culture 

; cultural  diffusion  process 
ask  families  [ 

if  any?  my -group  [ 

let  consensual - culture  compute  - consensual - culture -within- group 
let  diff erent - cultural - traits  [ ] 
foreach  n- values  cultural -dimensions  [?] 

[ if  (item  ? culture)  !=  (item  ? consensual - culture)  [ set  diff erent - cultural - traits  fput  ? different- 
cultural  - traits  ] ] 

if  not  empty?  diff erent - cultural - traits 
[ if  (random- float  1 < p - cult - dif fusion) 

[ let  trai t - to  - copy  one-of  dif f erent - cultural - traits 

set  culture  replace-item  trait  - to  - copy  culture  (item  trai t - to  - copy  consensual  - culture  ) ] ] 

] 

] 

; cultural  mutation  process 
ask  families  [ 

if  random-float  1 < p - cult -mutation 

[ let  trai t - to -mutate  one-of  n-values  cultural  - dimensions  [?] 
let  new- trait  random  cultural - trai ts 

set  culture  replace-item  trait  - to -mutate  culture  new-trait 

ask  my-group  [ set  culture  replace-item  trait  - to -mutate  culture  new-trait  ] 

] 


end 


5243  1 1 1 241 


243  1 1 1 241 


3523352 


Subprocesses.CULTURAL 

CONSENSUS 


ticks:  535 


52431 11241 


52431 11241 


Identity  comparison  after  652  ticks 


Subprocesses.  REPRODUCE 


Agents  lose  one  of  its  members  (a  labor  unit)  each  time  the  total  acquired 
energy  is  below  survival  threshold.  In  the  same  way,  every  30  ticks,  a new 
member  is  born,  and  it  will  be  alive  until  total  acquired  energy  is  below  survival 
threshold.  In  this  way,  we  have  implemented  a determinist  population  growth 
mechanism  opposed  to  stochastic  mortality. 


Subprocesses 


to  replace 

; Replacements  uses  a roulette  wheel  rule  according  to  families'  capacities 

; oldage  deaths 

let  oldage - deaths  [ ] 

ask  families  [ 

set  age  ( age  + 0.5) 

if  ( age  >=  max- age  ) [ set  oldage - deaths  (fput  who  oldage - deaths ) ] 

] 

; new  families 

let  parents  descendant -parents  length  oldage - deaths 
foreach  parents  [ new- descendant  ? ] 
foreach  oldage - deaths  [ ask  family  ? [ die  ] ] 

set  n-deaths -by-old-age  (length  oldage - deaths ) 

; starving  deaths 
let  starving  - deaths  [ ] 
ask  families  [ 

if else  ( energy  <=  0.01  ) 

[ set  periods -without  - energy  (periods -without  - energy  + 1)] 

[ set  periods -without  - energy  0 ] 

if  ( periods -without  - energy  > max-periods -without  - energy)  [ set  starving  - deaths  (fput  who  starving  - deaths ) ] 

] 

foreach  starving  - deaths  [ ask  family  ? [ die  ] ] 

set  n-deaths -by- insufficient  - energy  (length  starving  - deaths ) 
set  parents  descendant -parents  length  starving  - deaths 
foreach  parents  [ new- descendant  ? ] 

End 

to  new- descendant  [ parent -who  ] 
ask  family  parent -who 
[ hatch  1 [ 

set  age  0 

set  max-age  random -poisson  life  - expectancy 
set  capacity  random- float  1 
set -neighborhood -and -group  ] 


end 


] 


Experimental  Results 


Scenario  1 


Scenario  2 


Scenario  3 


Scenario  4 


Scenario  5 


Controlled  Variable: 

Population 


Controlled  Variable: 

Cooperation  Advantages 


Controlled  Variable: 

Cultural  Dimensions 


Controlled  Variable: 

Tollerance  to  others 
difference 


Controlled  Variable: 

Size  of  Mobility 
Neighbourhood  Area 


Scenario  1 


Experimento  1:  Media  de  Grupos  Formados 
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Experiment  1 
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Scenario  2 


Controlled  Variable: 

Size  of  Mobility 
Neighbourhood  Area 


Scenario  3 


Exp  2-  Variacion  de  lo  retornos  de  la  cooperacion 
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Retornos  de  la  cooperacion 
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Controlled  Variable: 

Cooperation  Advantages 


Experimento  2-  Variable  Retornos  de  la  cooperacion 
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Retorno  de  la  Coop  2 
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Groups  evolution 
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'he  consequences  of  cooperation:  diffusion  of  technological  innovations, 
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Scenario  4 


Relacion  entre  el  N°  de  grupos  formados,  dimensiones 
culturales  y n°  de  poblacion 


n°  de  agentes 


-2  Dimensiones  culturales 
6 Dimensiones  Culturales 
- 10  Dimensiones  Culturales 


— ■ — 4 Dimensiones  Culturales 
8 Dimensiones  Culturales 


Controlled  Variable: 

Identity 
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Scenario  5 


Relacion  entre  el  numero  de  poblacion  yel 
porcentaje  de  Tolerancia  Identitaria 
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Controlled  Variable: 

Intolerance  to  others 
difference 


efectos  en  la  variacion  de  los  retornos  de  cooperacion 


retorno  de  la  cooperacion 
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Scenario  5 


Groups  evolution 


0 time  20500 


Groups  distribution 

25 


Controlled  Variable: 

Intolerance  to  others 
difference 


Scenario  5 


Controlled  Variable: 

Intolerance  to  others 
difference 


The  consequences  of  cooperation:  social  aggregation. 


Max.  Resources  15000.  Max.  Resources  20000.  Max.  Resources  25000. 


Max.  Resources  40000. 


Max.  Resources  50000. 


The  consequences  of  cooperation:  social  aggregation 


Fractionalization, 

Herfindhal  Concentration  Index 


i'  = l /. 


Scenario  1 

Poor  resources 

Scenario  2 

Scenario  3 

Scenario  4 

Scenario  5 

Abundant 

resources 

ELF 

0.9863 

0.92 

0.97 

0.423 

0.431 

The  consequences  of  cooperation:  social  aggregation. 


Generalized  Resemblance 

Bossert  et  al.  Ethno-Linguistic  Fragmentation 


G(Sn)  = 1 


Scenario  1 

Poor  resources 

Scenario  2 

Scenario  3 

Scenario  4 

Scenario  5 

Abundant 

resources 

G(SN) 

4.96 

0.86 

0.06 

0.115 

0.144 

The  consequences  of  cooperation:  social  aggregation. 


Generalized  Resemblance 


Resources  15000 


Resources  20000 


Resources  25000 


Resources  40000 


Resources  50000 


The  consequences  of  cooperation:  social  aggregation 


Demographic  polarity 

Montalvo  and  Reynal-Querol  Index 
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Scenario  1 

Poor  Resources 

Scenario  2 

Scenario  3 

Scenario  4 

Scenario  5 

Abundant 

Resources 

RQ 

0.0511 

0.14 

0.09 

0.56 

0.56 

STATED  EFFECTS  OF  THE  SOCIAL  AGGREGATION 


High  benefits  of  cooperation,  since  the  households  will  cooperate  with  a 
major  quantity  of  households. 

Low  levels  of  tolerance.  Households  interact  less  with  those  households 
with  a minor  number  of  common  identity  features. 

Emerge  in  territorial  areas  with  higher  numbers  of  households  and/or  a higher 
density  of  population  provided  tolerance  levels  to  other  difference  remain 
relatively  low. 

The  higher  the  cultural  resemblance  and  proximity 

The  higher  the  size  of  the  neighborhood  and  the  size  of  the  population 


STATED  EFFECTS  OF  THE  SOCIAL  SEGREGATION 


•Less  benefits  of  labor  cooperation.  If  the  interaction  between  households  does 
not  generate  a progress  in  the  probabilities  of  survival,  households  will  survive 
isolated,  without  any  convergence  in  their  vectors  of  identity. 

•Low  density  of  population  and  low  intolerance,  given  the  low  probabilities  of 
social  interaction  and  labor  cooperation. 

•Even  in  case  tolerance  to  others  difference  be  low,  households  tend  to  fusion 
their  identities  when  the  number  of  shared  dimensions  overcomes  50  %. 

•Abundant  resources  and  a high  level  of  intolerance,  households  segregate 
easily. 

•When  the  size  of  neighborhood  area  increases,  and  the  level  of  intolerance 
overcomes  a threshold  of  50  %-60  %,  prior  aggregations  tend  to  split. 

•When  the  conditions  of  subsistence  deteriorate  (resources  diminish  and  the 
probabilities  for  survival  also  decrease)  prior  aggregations  tend  to  split. 


Everything  is  possible... 

if  you  just  j 

SIMM 

iTEj 

Not  everything  is  possible  in  a simulated  world 


What  is  possible  within  the  model? 

Sensitivity  Analysis 


Sensitivity  analysis:  Ensuring  that  the 
model  produces  reasonable  results 
when  factors  and  assumptions  are 
changed 

We  need  to  test  the  behavior  and 
response  of  the  model  when 
modifying  in  a directed  way  various 
assumptions  and  factors 


What  is  possible  within  the  model? 

Exploring  the  Parametric  Space 


In  sensitivity  analysis,  the  number  of  factors  should  be  kept  as  low  as  possible.  But  thi 


What  is  possible  within  the  model? 

Exploring  the  Parametric  Space 
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description 
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Model 

Implement 

ation 

(Netlogo) 
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RAW 

DATA 

What  is  possible  within  the  model? 

Exploring  the  Parametric  Space 


For  each  value  of  N (=10, 100,  200,  400), 
a repetition  has  3x4x3x3x3x4x3=  3888  runs 
(parameter  combination), 
each  one  stops  after  700  ticks. 


What  is  possible  within  the  model? 

Exploring  the  Parametric  Space 


Because  of  the  intrinsic  stochasticity  of  the  model 
executed  repeated  number  of  times  each  posible 


, we  have 
scenario 


We  have  made  30  runs  para  for  each  parametric  combination. 


Each  run  simulates  350  years  of  historical  evolutio 
each  cycle  represents  a season,  we  have  executed 
each  run. 


n.  Given  that 
700  cycles  for 
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possible 
within  the 
model? 


Exploring  the 
Parametric 
Space 
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Jj  N=10 
N=100 

N=400  (miercoles  24) 


N=10  -10  REP  [lunes  4 mayo) 
N=10  - 10  REP  [marten  5 mayo] 
N=10  - 10  REP  [vi ernes  1 mayo) 

N = 100  2 REP  [miercoles 20  mayo] 
N = 1003  REP  (jueves  14  mayo) 

N = 100  3 REP  [lunes  11  mayo] 

N = 100  3 REP  [vi ernes  22  mayo) 

N = 100  4 REP  [martes  12  mayo] 

N = 100  5 REP  [jueves  21  mayo) 

Nz  100 10  REP  (jueves  30  abril) 


What  is  possible  within  the  model? 

Exploring  the  Parametric  Space 

For  each  run  the  simulation  generates  4 files: 


ARCHIVOS  AGENTS 

ARCHIVOS  OUTPUT 

ARCHIVOS  SUPERAGENT 

PARAMETERS 

IUE: 

IUE 

IUE 

IUE 

cooperation - allowed? 

.ticks 

ticks 

ticks 

help -limited -to- spare - 

number- o f - f ami ly-de  a t hs - this-t i c k 

who 

labor? 

number-of-labor-units-deaths-this-tick 

initial-tick 

initial -population 

yjc or 

number rof -movements 

final -tick 

mean - resource - on - patches 

col 1 ec t ed-energy 

max- nurabe  r -o  f -membe  r s 

sd- resource-on-patches 

total -energy 

group-consensual-identity.  [01 

average -technology 

labor 

group-consensual-identity  [1] 

diversity 

survival - t hr esho Id 

qroup-consensual-identity  [2] 

labor -average 

technology 

average -storing -factor 

storing- factor 

group- consensual -identity  [ 3 1 

movement 

indi vi dua 1 - c apab i 1 i ty 

qroup-consensual-identity  [4] 

internal-change- rate 

s imi  1 a r i t y- thr  e s hold 

group -consensual- ident i t y [ 5 ] 

cooperation? 

group- consensual-ident i t y.  [ 6 ] 

tevrsuperagent 

qroup-consensual-identity  [7] 

Identity.  [01 

group  - con  sensual  - ident  i ty.  |[  8 1 

identity.  [1] 

qroup-consensual-identity.  [9] 

identity.  [2] 

identity.  [3] 

identity  [4] 

identity.  [51 

identity  [6] 

identity  [7] 

identity.  [8] 

Identity  [91 

What  is  possible  within  the  model? 

Exploring  the  Parametric  Space 

General  Structure  of  the  Agent  File:  each  record  is  a description  of  each  a 

IUE  ticks  who  xcor  ycor  collected-energy  total-energy  labor  survival-threshold  technology  storing-factor  individual-capability  similarity-threshold  cooperation?  my-superagent  identity  [0] 
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0.78 

0 false 

0 

4 

1E+08 

0 

5 

63 

20 

0 

1460.5 

4 

1460 

10.5 

0.648 

0 false 

0 

4 

1E+08 

0 

3 

28 

28 

0 

1095.5 

3 

1095 

10.5 
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8 

17 
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0 false 
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4 
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0 

8 

36 

13 
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8 
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0 false 
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6 

23 
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63 
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5 

1825 

10.5 

0.71 
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10 

4 
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1 

9 

56 

12 
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9 
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10 

4 

1E+08 
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6 

23 
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7 
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1 

46 

36 
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6 
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10 

4 
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1 

3 

28 

28 
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3 
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0,1  true 

10 

4 
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4 

43 

9 
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10 

4 
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0 
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39 
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10 
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36 

13 
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46 
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2190 
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56 

12 

3285 
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9 

3285 
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0.78 

0 false 

10 

4 
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4 

43 

9 
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2 
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0.424 

0 false 

10 

4 
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3 
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2 
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3 
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1E+08 
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0 
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2190 
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6 
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10.5 

0.634 

0 false 
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6 
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0 false 
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What  is  possible  within  the  model? 


Exploring  the  Parametric  Space 

General  Structure  of  the  Output  File:  Each  record  is  a description  of  an  individual  tick  fron 
of  living  agents  and  an  average  number  of  moving  agents  has  been  recorded. 


ticks 

number-of-family-deaths-this-tick 

number-of-labor-urits-deaths-this-tick 

number-of-movements 

100006656 

0 

0 

0 

0 

100006656 

1 

0 

0 

1 

100006656 

2 

0 

0 

0 

100006656 
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1 

100006656 

4 

0 

0 

1 

100006656 

5 

0 

0 

0 

100006656 

6 

0 

0 

0 

100006656 

7 

0 
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1 
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S 

0 

0 

0 

100006656 

9 

0 

0 

0 

100006656 

10 

0 

0 

0 

100006656 

11 

0 

0 

0 

100006656 

12 

0 

0 

1 

100006656 

13 

0 

0 

0 

100006656 

14 

0 

0 

0 

100006656 

15 

0 

0 

0 

100006656 

16 

0 

0 

0 

100006656 

17 

0 

0 

0 

100006656 

18 

0 

0 

2 

100006656 

19 

0 

0 

1 

in  n n n 


i nnnncccc 


What  is  possible  within  the  model? 

Exploring  the  Parametric  Space 


General  Structure  of  the  SuperAgent  File:  each  record  is  a descriptior 
a particular  moment  of  time  (tick)  for  each  run 


UE 

ticks 

who 

initial-tick 

final-tick 

max-number-of-members  group-consensual-identity  [0] 

group-consensual-identity  [1] 

group-consensual-identity  [2]  group-consensual-identr 

1E+08 
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10 

1 

0 

10 

4 

3 

1 

1E+08 

3 

10 

1 

2 

10 

4 

3 

1 

1E+08 

3 

11 

3 

0 

10 

4 

3 

1 

1E+08 

4 

11 

3 

4 

10 

4 

3 

1 

1E+08 

59] 
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59 

0 

11 

4 

3 
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IE +08 

60 

13 

59 

60 

11 

4 

3 

1 

IE +08 

871 

3 

87 

0 

12 

4 

3 

1 

1E+08 

88 

3 

87 

88 

12 

4 

3 

1 

1E+08 

131 

19 

131 

0 

15 

4 

3 

1 

1E+08 

132 

19 

131 

132 

15 

4 

3 

iT 

1E+08 

133 

20 

133 

0 

15 

4 

3 

iT 

1E+08 

134 

20 

133 

134 

15 

4 

3 

1 

1E+08 

137 

21 

137 

0 

15 
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3 

1 

1E+08 

138 

21 

137 

138 

15 

4 

3 

1 

What  is  possible  within  the  model? 

Big  Data  or  "a  lot  of"  Data 

We  assume  the  size  of  the  population  at  the  beginning  of  the  simulation  is  the  pan 
the  degree  and  nature  of  Cultural  Diversity  at  the  end  of  the  simulation 


Number  of  Generated  Files: 


1.394.205  files 


First  Scenario:  10  agents  at  start  up 

30  runs,  700  deles  simulated  at  each  run: 


116639  files  (Current  State  of  each  agent  at  each  tick  ) 


30  runs,  700  cides  simulated  at  each  run:  116639  files  (Number  of  living  Agents  by  tick) 

30  runs,  700  cicles  simulated  at  each  run:  96123  files  (Current  State  of  each  SuperAgent  at  each  tick 
30  runs,  700  cicles  simulated  at  each  run:  116639  files  (Global  parametres  of  each  run) 


Second  Scenario:  100  agents  at  start  up 

30  runs,  700  cicles  simulated  at  each  run: 


116639  files  (Current  State  of  each  agent  at  each  tick  ) 


30  runs,  700  cicles  simulated  at  each  run:  116639  files  (Number  of  living  Agents  by  tick) 

30  runs,  700  cicles  simulated  at  each  run:  116150  files  (Current  State  of  each  SuperAgent  at  each  tick 
30  runs,  700  cicles  simulated  at  each  run:  116639  files  (Global  parametres  of  each  run) 


Third  Scenario:  200  agents  at  start  up 

30  runs,  700  cicles  simulated  at  each  run: 


116639  files  (Current  State  of  each  agent  at  each  tick  ) 


30  runs,  700  cicles  simulated  at  each  run:  116639  files  (Number  of  living  Agents  by  tick) 

30  runs,  700  cicles  simulated  at  each  run:  116633  files  (Current  State  of  each  SuperAgent  at  each  tick 
30  runs,  700  cicles  simulated  at  each  run:  116639  files  (Global  parametres  of  each  run) 


Fourth  Scenario:  400  agents  at  start-up 

30  runs,  700  cicles  simulated  at  each  run:  Number  of  Agent  Files 

2 runs,  700  cicles  simulated  at  each  run:  3887  files  (Current  State  of  each  agent  at  each  tick  ) 

2 runs,  700  cicles  simulated  at  each  run:  3887  files  (Number  of  living  Agents  by  tick) 

2 runs,  700  cicles  simulated  at  each  run:  3887  iles  (Current  State  of  each  SuperAgent  at  each  tick 

2 runs,  700  cicles  simulated  at  each  run:  3887  files  (Global  parametres  of  each  run) 


A data  base  model  for  Millions  of  records 


Model 

Implement 

ation 

(Netlogo) 


ANALYSIS 

(R-studio) 


EXTRACTI 

ON 

(SQL 

language) 


Queries  in 
order  to 
make  data 
statistically 
manageable 


Model 

Run 
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RESULTS 


MySQL  Processing 

btaining  a DataBase  for  Statistical 

Analysis 


FOR  EACH  RUN  (IUE)  FOR  EACH  SELECTED  TICK  FROM  EACH  RUN 


IUE:  Number  of  Run 
TICK:  Selected  cicles  from 
each  run:  Select: 

0, 

1, 

50, 

51, 

100, 

101, 

200, 

201, 

300, 

301, 

400, 

401, 

500, 

501, 

698, 

699. 


GLOBAL  PARAMETERS  AT  START  UP 
COOPERATION  ALLOWED:  STATE  VARIABLE 
INITIAL  POPULATION  MEAN  OF  RESOURCES 
STAND.  DEVIATION  OF  RESOURCES 
SEASON:  Even  Number  of  Tick=  HOT; 

Odd  Number  of  Tick=  COLD 
INITIAL  TECHNOLOGY 
INITIAL  STORING  FACTOR 
DIVERSITY 

AVERAGE  LABOR  UNITS  PER  AGENT 
AVERAGE  POSSIBILITY  FOR  MOVEMENT 
RANDOM  CHANGE  RATIO  -ICR 


DEMOGRAPHY  AND  SURVIVAL 

NUMBER  OF  LIVING  AGENTS  AT  THE  CURRENT 
TICK 

TOTAL  NUMBER  OF  LABOR  UNITS  AT  THE 
CURRENT  TICK 

NUMBER  OF  LABOUR  UNITS  DEATHS  THIS  TICK 
AVERAGE  OF  COLLECTED  ENERGY  PER  AGENT 
AT  THE  CURRENT  TICK 

COOPERATION 

CURRENT  NUMBER  OF  COOPERATIVE  AGE 
CURRENT  NUMBER  OF  LIVING  SUPERAGEh 


GLOBAL  MEASURES  OF  CULTURAL  DIVERSITY 

SOCIAL  FRACTIONALIZATION 
DEPTH  OF  DIFFERENCES 
GENERALIZATION  OF  FRACTIONALIZATION 
GENERALIZED  RESSEMBLANCE 
GENERALIZED  SPATIAL  AGGREGATION 


AVERAGE  DURATION  OF  LIVING  AGENTS  A1 
CURRENT  TICK 

AVERAGE  SIMARITY-THRESHOLD  AT  THE 
CURRENT  TICK 


SOCIAL  INTERACTION 

MAXIMUM  VALUE  OF  TECHNOLOGY  BY  ANY  AGENT  AT  THE  CURRENT  TICK 
MINIMUM  VALUE  OF  TECHNOLOGY  BY  ANY  AGENT  AT  THE  CURRENT  TICK 
MAXIMUM  VALUE  OF  STORING  FACTOR  BY  ANY  AGENT  AT  THE  CURRENT  TICK 
MINIMUM  VALUE  OF  STORING  FACTOR  BY  ANY  AGENT  AT  THE  CURRENT  TICK 


MySQL  Processing 

3.  Obtaining  a DataBase  for  Statistical  Analysis 


Patagonial65  - Relational  DB  structure  and  keys 


parameters.txt 


IUE 


cooperation-allowed? 

he  I p-lim  ited -to-spare-l  abor? 

initial-population 

m ean-  reso  u rce-o  n -p  ate  h es 

average-technology 

diversity 

labor-average 

average-stori  n g-f  actor 

movement 

internal-change-rate 

displlay-resorces? 


Output.txt 


>IUE^ 


ticks 


nu  m be  r-of-f  am  ily- 
deaths-this-tick 

nu  m be  r-of -labor-u  ni  ts- 
deaths-this-tick 

num  be  r-of-move  merits 


agents.txt 


► IUE-* 


ticks  \+— 
who 
xcor 
ycor 

collected-energy 

total-energy 

labor 

survival -threshold 

technology 

storing -factor 

individual-capability 

similarity-threshold 

cooperation? 


my-superagent  ^ 


identity  (0-9) 


max-number-of- 

members 

group-consensual- 
identity  (0-9) 


MySQL  Processing 

3.  Obtaining  a DataBase  for  Statistical  Analysis 

FOR  EACH  SELECTED  TICK  FROM  EACH  RUN 


&QIBIC 


Limit  to  5H  rows 


I ^ 0,  HI  H 


agentsjtest.IUEj 
agents_test.ticfcsj 

Count ( age nts_test. who)  AS  CuentaDevihOj 
Sum ( age ntstest. labor)  AS  SumaDelaborj 
Avg (age ntstest. labor)  AS  Profiled ioDe labor 
FROM  agentsjtest 

GROUP  BY  agentstest.IUEj  agentsjtest. ticks 
ORDER  BY  agentstest.IUE*  agentstest. ticks; 


Result  Grid  | £§  ^ Hfter  Rows: 


0438411936 
0438411936 
0438411936 
0438411936 
0438411936 
0488411936  5 

0438411936  6 

0438411936  7 

0438411936  8 


MySQL  Processing 
3.  Obtaining  a DataBase 
for  Statistical  Analysis 


FOR  EACH  SELECTED  TICK  FROM  EACH  RUN 


SQL  QUERIES 


COOPERATION 

CURRENT  NUMBER  OF  COOPERATIVE  AGENTS 
CURRENT  NUMBER  OF  LIVING  SUPERAGENTS 

AVERAGE  DURATION  OF  LIVING  AGENTS  AT  THE 
CURRENT  TICK 

AVERAGE  SIMARITY-THRESHOLD  AT  THE 
CURRENT  TICK 


SELECT  DISTINCTROW 
superagents_test.iue, 
superagents_test.ticks, 
superagents_test.who 

FROM  superagents_test 
GROUP  BY  superagents_test.iue, 
superagents_test.ticks, 
superagents_test.who 

ORDER  BY  superagents_test.iue, 
superagents_test.  ticks; 


MySQL  Processing 

3.  Obtaining  a DataBase  for  Statistical  Analysis 


FOR  EACH  SELECTED  TICK  FROM  EACH  RUN  SQL  QUERIES 


COOPERATION 

CURRENT  NUMBER  OF 
COOPERATIVE  AGENTS 
CURRENT  NUMBER  OF  LIVING 
SUPERAGENTS 
AVERAGE  DURATION  OF 
LIVING  AGENTS  AT  THE 
CURRENT  TICK 
AVERAGE  SIMARITY- 
THRESHOLD  AT  THE  CURRENT 
TICK 


1 

IUE  ticks 

who 

initial-tick 

final-tick 

max-number-of-members 

2 

100006656  1 

10 

1 

0 

10 

3 

100006656  2 

10 

1 

2 

10 

4 

100006656  3 

11 

3 

0 

10 

5 

100006656  4 

11 

3 

4 

10 

6 

100006656  59 

13 

59 

0 

y 

Superagent  13  is 

7 

100006656  60 

13 

59 

60 

created  at  tick  59  and 

8 

100006656  87 

15 

87 

0 

it  dies  at  tick  60 

9 

100006656  88 

15 

87 

88 

12 

10 

100006656  131 

19 

131 

0 

15 

11 

100006656  132 

19 

131 

132 

15 

12 

100006656  133 

20 

133 

0 

15 

13 

100006656  134 

20 

133 

134 

15 

14 

100006656  137 

21 

137 

0 

15 

IS 

100006656  138 

21 

137 

138 

15 

16 

100006656  139 

22 

139 

0 

15 

17 

100006656  140 

22 

139 

140 

15 

18 

100006656  159 

25 

159 

0 

15 

19 

100006656  160 

25 

159 

160 

15 

20 

100006656  161 

26 

161 

0 

16 

21 

100006656  162 

26 

161 

162 

16 

22 

100006656  163 

27 

163 

0 

15 

23 

100006656  164 

27 

163 

164 

15 

24 

100006656  165 

28 

165 

0 

15 

25 

100006656  166 

28 

165 

0 

15 

26 

100006656  167 

28 

165 

0 

Y 

Superagente  28  is 

27 

100006656  168 

28 

165 

168 

created  at  tick  165  and 

28 

100006656  185 

32 

185 

0 

it  dies  at  tick  168 

29 

100006656  186 

32 

185 

186 

18 

30 

100006656  189 

33 

189 

0 

16 

31 

100006656  190 

33 

189 

190 

16 

MySQL  Processing 

3.  Obtaining  a DataBase  for  Statistical  Analysis 


FOR  EACH  SELECTED  TICK  FROM  EACH  RUN  SQL  QUERIES 

CREATE 

ALGORITHM  = UNDEFINED 
DEFINER  = 'newuser'@'%' 

SQL  SECURITY  DEFINER 
VIEW  'query  2 moviment  poblacio'  AS 
SELECT 

'agentsJest'.'IUE'  AS'IUE', 
'agents_test'. 'ticks'  AS  'ticks', 
MIN('agents_test'.'xcor')  AS  'xmin', 
MAX('agents_test'.'xcor')  AS  'xmax', 
MIN('agents_test'.'ycor')  AS  'ymin', 
MAX('agents_test'.'ycor')  AS  'ymax' 
FROM 

'agents_test' 

GROUP  BY 'agentsJest'.'IUE'  , 
'agents  Jest'. 'ticks' 

ORDER  BY  'agentsJest'.'IUE'  , 

'agents  Jest'. 'ticks' 


COOPERATION 

CURRENT  NUMBER  OF 
COOPERATIVE  AGENTS 
CURRENT  NUMBER  OF  LIVING 
SUPERAGENTS 
AVERAGE  DURATION  OF 
LIVING  AGENTS  AT  THE 
CURRENT  TICK 
AVERAGE  SIMARITY- 
THRESHOLD  AT  THE  CURRENT 
TICK 


MySQL  Processing 
3.  Obtaining  a DataBase 
for  Statistical  Analysis 


FOR  EACH  SELECTED  TICK  FROM  EACH  RUN 


SOL  GUERIES 


rar  - View  query  4 - agents  i labor  SQ  L File  24* 

S Q uery7-tech  nology-  i-stori  ng  - V...  Q uery7-tech  nology-  i-stori  ng 

\a  b \t  f it o \m  e 

Limit  to  50  raws  T | -0  fij]  |?| 

SOCIAL  INTERACTION 

MAXIMUM  VALUE  OF  TECHNOLOGY  BY  ANYAGEN 
MINIMUM  VALUE  OF  TECHNOLOGY  BY  ANY  AGEN 
MAXIMUM  VALUE  OF  STORING  FACTOR  BY  ANYA 
MINIMUM  VALUE  OF  STORING  FACTOR  BY  ANY  AG 


1 & 

2 

3 

4 

5 

6 
7 
B 
9 

10 

11 


SELECT 

'agents  test' .'IUE'  AS  IUE', 

' agent stest' .'ticks'  AS  'ticks', 

MIN (' agent s_test' .'technology' } AS  'MinDetechnology' , 
'technology')  AS  'MaxDetechnology' , 
'storing')  AS  ' Mi n&e storing' , 
'storing')  AS  'MaxDe storing' 


MAX(' agentsjtest' 

MIN(' agents_test' 

MAX(' agents_test' 

FROM 

' agent s_test' 

GROUP  BY  ' agentsjtest' .'IKE' 
ORDER  BY  ' agentsjtest' .' IUE' 


' agentsjtest' . 'ticks' 
' agentsjtest' . ' ticks'l 


l|  ResuItGrid  _[|j 

Fter  Rows: 

Export:  Wrap  Cel 1 

Content:  Fetch  rows:  EH?  ^ 

□ 

IUE 

ticks 

MinDetechnology 

MaxDetechnology  MinDestoring 

MaxDestoring 
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FI 

0438411936 
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0,01 
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043S411936 
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0.582 

2 

0,01 

0.91 

0438411936 

3 

0.582 

2 

0,01 

0.91 

0438411936 

4 

0.582 

2 

0,01 

0.91 

0438411936 

5 

0.582 

2 

0,01 

0.91 

0438411936 

6 

0.588 

2 

0,01 

0.91 

0438411936 

7 

0.588 

2 

0,01 

0.91 

0438411936 

8 

0.588 

2 

0,01 

0.91 

MySQL 

Processing 


1.  File 

Aggregation 
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Four  files 

Each  record 

is  a description  of  each  agent 
at  each  tick  (Cycle)  for  each  run. 


Initial  population  =10 
3 million  RECORDS 


Initial  population  =100 
More  than  3 million  RECOF 


Initial  population  =200 
More  than  5 million  RECOF 


Initial  population 
=400 

700  records 


B 


| IUE  ^ IUE  ^ IUE 
0438411336  0438414579  0438425497 


ticks 


10 


11 


12 


13 


14 


IS 


19 


20 


21 


22 


23 


24 


25 


26 


27 


who 


who 


13 


14 


14 


14 


14 


14 


14 


14 


14 


14 


14 


14 


14 


14 


14 


14 


14 


14 


14 


14 


14 


14 


14 


14 


14 


14 


14 


14 


who 


10 


10 


10 


10 


IUE 


0438430482 


who 


10 


700  ticks 


350 


-0438411936 


-0438414579 


-0438425497 


-0438430482 


97  ticks 


20 


0438411936 


1 0438411936  1 0438414579  1 0438425497  1 0438430482 


We  have  selected  c 
Cycle  (warm  seasoi 
one  cycle  (cold  sea 
every  100  cycles. 


Resulting  file  has 
300000  records. 
For  each  N=  10, 
N=  100, 

N=  200, 

N=  400 


A data  base  model  for  Millions  of  records 


Model 

Implement 

ation 

(Netlogo) 


ANALYSIS 

(R-studio) 


EXTRACTI 

ON 

(SQL 

language) 


Queries  in 
order  to 
make  data 
statistically 
manageable 


Model 

Run 

(Netlogo) 


DDBB 

(MYSQL) 


RAW 

DATA 


i 


Hundreds  of 
thousands 
files 
Approx. 
500.000 


INTEGRAT 

ION 


Classify 

data 


Hundreds  of 
millions  records 
Approx.. 
400.000.000 


RESULTS 


Value 


Exploring  the  Parametric  Space 

The  sources  of  Stochasticity 


Before  answering  why  each  run  gives 
Different  results,  we  have  to  measure 
the  global  stocasticity  of  the  system:,  and  to 
Ask  ourselves  about  the  consequences  of: 

Agents  find  other  agents 
just  by  chance. 

Individual  agents  have  the  the 
possibility  of  refuse  cooperation 
without  “rational”  reasons 
for  doing  it. 

Agests  “evolve”  also  for  particular 
reasons  and  not  only  as  a 
consequence  of  cooperation. 


Exploring  the  Parametric  Space 

The  Effects  of  Population  size  at  start-up 


Expected  Result:  the 
higher  the  population  at 
start-up,  the  “easier”  to 
build  cooperation 
networks,  the  more 
cooperation,  and  the 


Nevertheless,  the  higher  the  population  at 
he  population  will  grow  more,  and  control 
will  become  more  relevant. 
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Exploring  the  Parametric 

Space 

The  Effects  of  Resource  Availability 

The  Rich  World  / Poor  World  Scenarios 


Expected  Results:  The  poorer  the 
world,  the  higher  the  cooperation 
benefits,  and  the  more  stable  will  be 
cooperation  networks 


Nevertheless,  the  poorer  the  world,  the  higher  the 
Number  of  deaths,  and  the  less  population  density, 
what  is  a handicap  for  complex  cooperation  networks 


Why  Hunter-Gatheres  did  not 


Exploring  the  Parametric  Space 

The  advantages  of  cooperation  in  mobile  hunter-gatherers 


The  similarity-threshold  take  values 
between  0 and  1 , and  came  to  represent 
how  similar  should  be  an  agent  to  receive 
some  help.  It  simulates  the  “rationality”  in 
the  decision  whether  to  cooperate. 

My  similarity-threshold  depends  on  many 
factors  (number  of  labors,  level  of 
resources  in  my  patch,  difficulty  to  get 
them  ...)  but,  in  short,  is  proportional  to 
how  difficult  it  will  be  for  me  to  reach  my 
threshold  of  survival  in  the  current  tick. 

There  are  two  extreme  cases:  if  the  agent 
is  going  to  achieve  survival,  my  similarity- 
threshold  is  0 (always  help),  and  if  the 
agent  believes  its’s  not  going  to  reach 
(need  to  ask  for  help)  my  similarity 
threshold  is  1 . 
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Exploring  the  Parametric  Space 

The  advantages  of  cooperation  in  mobile  hunter-gatherers 
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The  decisive  question  in  any  simulation  is 
what  evidence  one  has  to  judge  the 
simulation  to  be  true. 


■'You  are  completely  tree  to  carry  out  whatever  research 
you  want,  so  long  as  you  come  to  these  conclusions/1 


Sensitivy  Analysis  is  not  Validation:  We  need  comparing 
simulated  results  with  empirical  data 


COMPARISON 
OF  COMPUTATIONAL 
MODELS 
O DOCKING 


THEORETICAL 

ANALYSIS 


LINKS  WITH 
THE  REAL 
WORLD 


CALIBRATING  THE  GLOBAL  THEORETICAL  MODEL  WITH 
PARTICULAR  ETHNOHISTORICAL  INFORMATION 


Number  of  Families 

[40] 

Size  of  Neighborhood  Area 

[5] 

Initial  Area  of  Settlement 

[50] 

Intolerance  Level 

[0-2] 

Num.  of  Cultural  Dimensions 

[10] 

Num.  of  Cultural  Features 

[10] 

Probability  of  cultural  diffusion 

[0.7] 

Probability  of  cultural  change 

[0.5] 

Cooperation  Advantages 

[4] 

Subsistence 

[0-2] 

Energy  Depretiation 

[1.7] 

Life  Expectancy 

[70] 

COMPARING  SIMULATED  RESULTS  WITH  ETHNONYMS 

(INDIGENOUS  POPULATION  SELF-IDENTITY) 


COMPARING  SIMULATED  RESULTS  WITH  ETHNONYMS 

(IMPOSED  ETHNICITIES) 


COMPARING  SIMULATED  RESULTS  WITH  LINGUISTIC 

DISTANCES 


COMPARING  SIMULATED  RESULTS  WITH 

ARCHAEOLOGICAL  DATA 

1.  Genetic  Distances  between  populations 
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COMPARING  SIMULATED  RESULTS  WITH 

ARCHAEOLOGICAL  DATA 
2.  Archaeological  Similarity  between  lithic  tools 


COMPARING  SIMULATED  RESULTS  WITH 

ARCHAEOLOGICAL  DATA 
3.  Archaeological  Similarity  between  Rock  Art  sites 


COMPARING  SIMULATED  RESULTS  WITH 
ARCHAEOLOGICAL  DATA 

3.  Prehistoric  Exchange  Networks:  the  source  of  raw  materials 


Patagonia  Simple  Past 


Limitations  and 
Enhancement 


Scale  of  Analysis:  Population  dynamics 

Z>  Scale  of  Analysis:  Temporal  units 

2>  Input  variables:  “Cultural”  Dimensions. 
Measuring  “Identity” 

Output  variables:  Number  of  Culturally  Similar 
Groups. 

How  to  measure  cultural  similarity.  Index  of  Ethnic 
fragmentation.  Degree  of  Social  Inertia  between  states 


5)  PROBLEMS  WITH  AGENT  RATIONALITY:  optimal 
choice,  bounded  rationality  or  bayesian  probability. 
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"Hunting,  gathering....  It's 
so  hard  to  prioritize!" 


6)  ENHANCING  THE  MODEL,  IMPLEMENTING  NEW 
MECHANISMS:  Other  mechanisms  of  cooperation 
restriction 

a.  Landscape  impact  on  mobility: 

i.  Atraction  forces:  rich  resources 

ii.  Repulsion  forces:  poor  resources  (negative).  Topographic  barriers  are  areas 
with  negative  resources 

b.  Social  impact  on  mobility:  kinship  vs.  Neighborhoods 

i.  Attraction:  alliances 

iii.  Repulsion:  enemies 

c.  Political  impact  on  mobility: 

i.  Attraction  forces:  leadership 

iv.  Repulsion  forces:  carrying  capacity,  social  conflict 

d.  Identity  similarity  vs.  Biological  phenotipic  similarity: 
from  ethnicity  to  race 


Recreating  the  Past  in  the 

Computer 
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“Would  someone  please  remind  me  what  our  origina 

intention  was?" 


Replication 


Simulation 

results 


In  archaeology,  agent-based 
simulations  have  been 
used  to 

• to  test  the  internal 
coherence  of  social  theories 
of  human  action  in  the  past 

• Explaining 
historical/archaeological 
data  with  some  hypothesis 
from  ethnographic  research. 


Simulating  animal  behavior 


Population  Number 


Animal  behavior  is  a good  example  of  social  mechanism  (without 
abstract  beliefs  nor  complex  motivations,  nor  desires  and  only 
simple  instinctive  intentions),  and  therefore  it  has  been  studied  in 
formal  terms  since  the  times  of  A.J.  Lotka  (1910)  and  V.  Volterra 
(1926). 


FROM  ANIMALITY  TO 
HUMANITY. 


There  is  a lot  of  “animality”  within 
us,  and  if  we  want  to  know  why 
we  do  what  we  are  doing  in  the 
present,  the  only  way  is  to 
understand  our  “degree  of 
animality”  and  the  historical 
process  of  differentiation  from 
our  “original”  animality. 


In  any  case,  the  most  important 
aspect  of  investigation  will  not  be 
the  animal  basis  of  human 
behavior,  but  the  specific 
process  of  progressive 
differentiation  in  the  way  we  take 
decisions  -more  or  less  rational- 
from  the  original  animal  instincts. 


FROM  ANIMALITY  TO 
HUMANITY. 


Among  recent  essays  in  this  direction,  we  can  mention:  Arenas  (2012), 
Hoban  et  al.  (2012),  Kawecki  (2012),  Ma  et  al.  (2012),  Kutsukake  and  Innan 
(2013),  Messer  (2013),  Mode  et  al.  (2013),  Villmoare  (2013),  Schlotterer  et 
al.  (2014),  Smaldino  et  al.  (2013),  Acevedo-Rocha  et  al.  (2014),  Hunemann 
(2014),  Lehman  and  Stanley  (2014),  Vevgari  and  Fioley  (2014),  Roseman  et 
al.  (2015),  Peart  (2015),  Shamrani  et  al.  (2015),  Smith  et  al.  (2015). 


FROM  ANIMALITY  TO 
HUMANITY. 
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WHAT'S  IT 
ALL  ABOUT? 
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In  any  case,  natural  selection  and 
evolutionary  mechanisms  have 
affected  animals  and  humans  not 
only  in  morphology,  but  in  the 
development  of  pre-human  behavior. 
It  is  also  the  question  of  the  origins 
of  “intelligence”  and  complex 
decision  making  (Gabora  and 
Russon  2011,  Gabora  and  DiPaola 
2012,  Kurzweil  2012, 
Chandrasekaran  2013,  Pringle  2013, 
Guddemi  2014,  Ross  and  Richerson 
2014,  Geary  2015)  and  also  culture. 
This  is  not  the  place  to  define  what  is 
culture,  but  recent  work  suggests  its 
computable  basis 


Out  of  Africa  Hypothesis 


The  causes,  mechanism  and 
Routes  of  expansion  of  Homo 
Erectus  have  been  analysed 


OTHER  EXAMPLES 


A computer  model  of  wayfinding 
behaviour  for  studying  hominin 
dispersals  and  their  chances  of 
longterm 
success. 

Added  complexity  comes  from 
considering  an  agent's  wayfinding 
motivation,  strategy,  and  differential 
cognitive  spatial  abilities 


Neanderthal  Behavior 


It  can  be  of  interest  to  compare  the 
dispersal  mechanism  of  pre-humans, 
to  the  motivations  and  intentionality 
of  movement  and  dispersal  by 
modern  humans  of  “prehistoric” 
times,  with  motivations  different  from 
modern  humans  of  present  times, 
and  even  our  antecessors  from  a 
more  recent  past  with  motivations 
assumed  to  be  like  ours 


Hunter-Gatherers 


GHA<  WAS  A HUNTER,  ZH00M8A  WAS  A GATHERER 
ANP  FNOK  WAS  A SPONGER. 


HUNTING-AND-GATHERING  IN 
THE  PAST  EXPLAINS  HOW  WE 
HAVE  SURVIVED  UNTIL  THE 
PRESENT. 

Prehistoric  hunter-gatherers  have 
been  studied  many  times  from  the 
point  of  view  of  animal  foraging 
behavior,  stating  that  human  agents 
also  forage  in  such  a way  as  to 
maximize  their  net  energy  intake  per 
unit  time. 


Hunter-Gatherers 
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The  understanding  of  many 
ecological  concepts  such  as 
adaptation,  energy  flow  and 
competition  hinges  on  the  ability  to 
comprehend  what  food  items  such 
human  agents  selected,  and  why. 
Nevertheless,  it  is  obvious  that  if 
humans  were  in  the  past  just  like  any 
other  animal  forager  or  predator,  we 
would  say  that  prehistoric  hunter- 
gatherers  survival  would  have 
depended  just  on  the  availability  of 
edible  resources.  Given  what  we 
know  about  the  natural  irregularity  of 
natural  resources  yield,  Homo 
sapiens  would  have  extinguished 
many  times  since  their  African 
origins! 


Hunter-Gatherers 


The  second  part  of  the  problem  when 
trying  to  couple  the  social  and  the 
environmental  lies  in  modeling  carrying 
capacity  and  the  capability  of  prehistoric 
humans,  even  with  inefficient  technology 
to  alter  and  modify  it. 

In  the  process  of  truly  coupling  human 
activity  and  natural  environment, 
computer  simulation  approaches  allow 
understanding  how  human  decisions 
and  subsequent  actions  would  change 
(at  least  affect)  the  structure  and 
function  of  many  natural  systems.  Such 
structural  and  functional  changes  would 
in  turn  exert  influence  on  human 
decisions  and  actions 


Reformulating  Malthus 


Old  Malthusian  view  on  population 
increasing  exponentially  while  food 
production  would  have  increased 
only  linearly,  in  constant  increments 
has  been  reformulated 
within  ABM  models  (Portugali  1999, 
Read  and  LeBlanc  2003,  Lane  2010, 
Cai  2012,  Schlueter  et  al.  2012, 

Levin  et  al.  2013,  Hritonenko  and 
Yatsenko  2013,  Ribeiro  2015). 


Time 


Prehistoric  Demography 


"We  have  one  hunter  and  or>e  gatherer...  everyone 
else  is  a consultant." 


The  single  most  obvious  constraint  of  human 
action  in  a particular  environment  is 
population  size,  especially  when  the  means  of 
production  seem  to  be  underdeveloped 
(hunting-and-gathering).  Many  modern 
computer  simulations  on  human  demography 
are  centered  on  modeling  the  particular 
dependence  on  annual  fertility  tables  and 
adopt  a fecundity  based  model.  The  odds  of 
conception  for  any  one  mating  event  can  be 
kept  constant  for  a female  agent  of  a given 
age,  and  the  probability  of  reproduction 
therefore  becomes  dependent  on  the 
frequency  and  timing  of  the  female  agent’s 
mating  activity.  This  allows  for  realistic  fertility 
variations  as  a function  of  mating  behavior 
frequency  (and  thus  contextual  opportunity  in 
the  form  of  access  to  male  sexual  resources) 
and  the  variations  of  individual  agent 
fecundity  over  time. 


The  myth  of  the  stupid  prehistoric  savages. 


Socio-ecological  models  make  emphasis  on 
physiological  motivation.  Intelligence  is 
expressed  in  the  way  they  look  for  the 
satisfaction  of  their  full  stomachs.  However,  if 
physiological  motivation  is  the  only  source  of 
directness  in  the  computer  simulation  of 
human  behavior  we  may  end  with  undesired, 
uniform  behavior. 

Human  (and  even  animal)  rationality  is  much 
more  complex  than  expected  and  therefore,  it 
is  easy  to  conclude  that  deterministic 
relationships  between  environmental  stress 
and  social  change  are  inadequate 


The  myth  of  the  stupid  prehistoric  savages. 


Socially  intelligent  agents  (SIAs)  should  be 
defined  as  agents  that  do  not  only  from  an 
observer  point  of  view  behave  socially  but 
that  are  able  to  recognize  and  identify  other 
agents  and  establish  and  maintain 
relationships  to  other  agents  (Dautenhahn 
1998).  The  process  of  building  SIAs  will 
always  been  influenced  by  what  the  human 
as  the  designer  considers  “social,”  and 
conversely,  agent  tools  that  are  behaving 
socially  can  influence  human  conceptions  of 
sociality 


The  myth  of  the  stupid  prehistoric  savages. 


One  of  the  very  first  computer  simulations  of 
prehistoric  hunter  gatherers  was  that  of 
Robert  Reynolds  (1986).  He  explicitly 
approached  the  problem  of  rationality  in 
hunter-gatherer  decision-making  in  terms  of: 
the  ability  of  each  member  to  collect  and 
process  information  about  the  resource 
distribution, 

the  extent  to  which  information  is  shared 
among  members, 

the  specific  sets  of  decision  available  to  each 
member,  and 

the  way  in  which  the  individual  decisions  are 
integrated  to  produce  a group  decision. 


Jim  Doran 
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: EOS  Model 


1'HE  EOS  MODKI 


local  density  of 
> agent  population 


individual  agent's 
subsistenoe  behaviour 


Key 


organization  and  co-operative 
problem  solving  effectiveness 
of  agent  population 


■npacts  upon 


Figure  9.2  I he  heart  of  the  KOK  model.  Population  concentration  impacts  upon 
the  characteristics  of  the  environmental  resource  acquisition  task,  and  both  impact 
on  agents*  social  models.  These  in  turn  affect  the  agents’  behaviour  and  hence  com- 
munity effectiveness,  which  lias  an  effect  on  on  population  density  closing  the  loop. 


THE  MYTH  OF  THE  GOOD 
PREHISTORIC  SAVAGE: 


The  origins  of  social 

differentiation  and  complexity. 

An  obvious  consequence  of  this 
way  of  considering  social 
organization  as  an  emergent 
property  of  the  mechanisms  of 
cooperation  (or  the  lack  of  it)  and 
cultural  transmission,  is  that  the 
origins  of  social  diversity, 
hierarchy  and  complexity  can 
also  be  considered  as  emergent 
properties  of  relatively  basic 
social  mechanisms. 


SIMULATING  SOCIAL 
TRANSITIONS 


The  Origin  of  Agriculture 

This  is  the  obvious  domain  for  computer 
simulations.  Can  agriculture  and  related 
practices  of  animal  control  emerge 
“mechanically”  in  a group  of  agents 
originally  defined  as  foragers  and 
predators?  The  technological  side  of  this 
transition  is  not  the  result  of  the 
“intelligence”  of  some  individuals  who 
“invented”  something  new.  As 
computational  simulations  have  proved 
domestication  of  plants  and  animals  is  an 
evolutionary  emergent  result.  Therefore  it 
seems  that  there  are  some  possibilities 
that  one  of  the  most  relevant  transitions 
in  the  history  of  humanity  had  also  a 
mechanical  basis. 
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Once  human  communities  begin 
producing  food  and  exchanging 
goods,  social  dynamics  become  more 
complex. 

There  are  a lot  of  examples  of  ABM 
models  of  complex  prehistorical 
societies 
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Different  authors  have  generated  computer  simulations  to 
understand  how  hierarchical  decision-making  could  have 
affected  inter-group  conflicts  sometime  through  the 
historical  evolution  of  human  society  (Mark  1998, 
Suleiman  and  Fischer  2000),  the  dynamics  of  status 
symbols  in  hierarchically  ordered  societies  (Pedone  and 
Conte  (2001),  the  consequences  of  wealth  distribution 
(Impullitti  and  Rebmann  2002),  the  coevolution  of 
farming  and  private  property  (Bowles  and  Choi  2013, 
Cockburn  et  al.  2013,  Angourakis  et  al.  2014),  the  origins 
of  war  (Duering  and  Wahl  2014)  and  the  Neolithic 
transition  from  egalitarianism  to  leadership  and 
despotism  (Levine  and  Modica  2013,  Powers  and 
Lehman  2015). 


OTHER  EXAMPLES 


A computer  simulation  of  the  Bronze 
Age  salt  mines  of  Hallstatt  in  Austria 
(dated  to  1458-1245  B.C.) 

By  employing  certain  hypotheses 
arising  from  excavated  archaeological 
findings  (e.g.  pick  handles,  lighting 
chips,  mine 

timber)  together  with  additional 
assumptions  and  estimated  data  (e.g. 
crop,  spatial  distances,  working  time), 
we  can  develop  a 

technological  reconstruction  in  form  of 
a mental  model. 
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OTHER  EXAMPLES 
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A computer  implementation  of 
mechanisms  of  population  dynamics 
and  carrying  capacity  of  the 
late  Iron  Age  oppida  in  Bohemia. 

The  system  dynamics  component  is 
used  to  provide  a way  to  explore  and 
test  various  general  theoretical 
hypotheses  related  to  the  functioning 
of  the  settlements  and  the  societal 
rules  that 
are  shaping  them, 


OTHER  EXAMPLES 
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A model  of  the  the  spatial  distribution 
of  early  Etruscan  states  within  a 
geographic  network  structure  in  order 
to  shed  light  on  how  space 
conditioned  their  development. 

The  relative  spatial  interactions 
amongst  the  sites  in  the  study  area 
has  been  simulated  using  an 
agentbased  model,  a radiation 
analogy  model,  and  a Hamiltonian 
gravitational  model. 


OTHER  EXAMPLES 
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A computational  model  of  hypothetical  trade  mechanisms  in  the  Roman 
East. 

The  Roman  Bazaar  is  based  on  the  idea  that  markets  in  the  Roman 
Empire  were  very  weakly  integrated  and  largely  driven  by  the  particular 
commercial  opportunities  of  agents  in  the  many  trading  places  around  the 
Empire.  These  opportunities  were  structured  by  the  agents'  social 
networks,  leading  to  a relative  unpredictability  of  supply  and  demand,  high 
uncertainty  of  information,  and  many  irregularities  in  prices  and  available 
products. 


OTHER  EXAMPLES 


A computer  model  of  the  settlement 
process  of  Oceania  based  upon  the 
simulation  of  trajectories  of  different 
types  of  canoes  sailing  through 
different  wind  patterns  throughout 
Oceania  to  suggest  that  there  were 
two  “schools  of  navigation”,  one  to  the 
east  and  another  to  the  west  of 
Samoa. 


OTHER  EXAMPLES 


Gavrilets  et  al.  (2014)  have 
developed  a spatially  explicit  agent- 
based  theoretical  model  of  the 
emergence  of  early  complex  polities 
via  warfare.  In  this  model  polities  are 
represented  as  hierarchically 
structured  networks  of  villages  whose 
size,  power,  and  complexity  change 
as  a result  of  conquest,  secession, 
internal  reorganization  (via  promotion 
and  linearization),  and  resource 
dynamics. 


OTHER  EXAMPLES 


A possible  criticism  about  the  idea  of 
simulating  the  past  would  be  the 
impossibility  of  simulating  the  historical 
evolution  of  complex  polities  in  modern 
times  for  reasons  of  scale:  to  be  fully 
capable  of  understanding  historical 
dynamics  of  ancient  empires  and  modern 
nations  we  would  need  to  create  artificial 
societies  of  such  complexity  that  any 
computer  could  run  the  simulation. 
Nevertheless,  the  current  use  of  agent-base 
modeling  and  related  techniques  to 
understand  modern  economics  and  modern 
social  and  political  organization  clearly 
indicates  the  opposite.. 
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Possible  Criticisms 


Are  we  LEARNING  something  about 
the  past 

or  just  playing  with  the  results  of 
our  historical  imagination? 


Possible  Criticisms 


We  have  to  investigate 
The  proper  relationship  between 
Historical  hypothesis 
Social  mechanism 
Archaeological  data 
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In  that  sense, 

Asocial  mechanism  without 
Temporal  dynamics  has  not  meaning 
For  a scientific  archaeology 
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But  Temporal  dynamics 
Can  only  be  defined  on 
The  basis  of  well  dated 
Historical-archaeological 
data 


Remember 


we  may  study  a subset  of  social  activity: 

collective  action, 


The  negative  side  of  this  approach  is  that 
there  is  no  possibility  of  knowing  why  an 
individual  person  made  something  somewhere 
at  some  moment.  However,  it  does  not 
presuppose  the  implicit  randomness, 
subjectivity,  or  indeterminism  of  social  action. 

The  goal  should  be  to  explain  the  sources  or 
causes  of  that  variability,  and  not  exactly  the 
inner  intentions  of  individual  action. 


Social  mechanisms 


It  is  not  complexity  in 
itself  what  should 
interest  us. 

Asocial  mechanism  is 
not  “better”  than  other 
because  it  is  more 
complex, 

But  it  is  its  explanantory 

capability 
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EXPLANANTION 


...  and 

“explanatory”  is 
not 

a universal 
property  of  any 
system.  It  is 
always  relative 
to  concrete 
circumstances, 
to  particular 
data  sets 
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Prediction 


"r  remember  when  tea  leaves  were  the 

STATE-OF-THE-ART  METHOD  OF  FORECASTING/' 


Computer  simulation  of  the  past 
Can  be  seen  as  a reivindication  for 
historical  prediction. 

If  we  feed  a computer  with  data 
from  the  past,  can  we  predict  the 
future? 
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The  answer  is 
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Oh,  not  a computer,  B-ctrcel 6 — 
we're  replacing  you  with  a 
computer  simulation. " 


